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THEORY OF THE ORIGIN OF THE VERY LOW-FREQUENCY RADIO EMISSIONS 
FROM THE EARTH’S EXOSPHERE* 


J. W. MacArthur 
Marlboro College, Marlboro, Vermont 
(Received May 13, 1959) 


In the theory propounded by Gallet’ to explain 
dawn chorus, the electromagnetic waves are 
generated by interaction between the natural 
thermal noise present in the exosphere and the 
beams of high-speed protons incident upon the 
earth’s magnetic field from the sun. This in- 
teraction, of the travelling-wave tube type, re- 
quires that the electromagnetic wave and the 
moving ions have essentially the same velocity 
over an extensive region of space in order for an 
appreciable amount of energy to be transferred 
from the beam to the wave. The following theory 
removes this objection, yet gives the same re- 
sults. 

A cloud of protons, impinging on the earth’s 
magnetic field at high speed, will in general be 
deflected into a helical path about the lines of 
force, rotating with angular frequency w, =He/Mc, 
where the proton mass M=1836 m,. Thus a 
proton gyrofrequency is set up which is 1/1836 of 
the electron gyrofrequency 

Wy =He /m_,c. 

This rotating cloud of protons may be consid- 
ered as a source of electromagnetic waves having 
angular source frequency w, and moving forward 
vith velocity V. Hence a Doppler shift in fre- 
quency will result, with the radiated angular fre- 
quency given by 


w=w,/(1-V/v,), (1) 
where Up» the phase velocity of the wave, is 
% =c/n. (2) 


Now, the index of refraction, m, in the “whistler 
mode” of propagation through the exosphere, is 
known’ to be given by 


W=l+w*/(ww, -w*), (3) 
where wy = (4nNe*/m,)”*, the plasma angular 


frequency. 
Combining (1), (2), and (3), there results 


w*(1 - V*/c?) - w*[w, - (V*/c*) Ww, +2. | 
+ w[w.* +2w w+ (V2 /c*)w,”] - (ww) =0. 
If V*<< c? and ws << wy, which are both very 
good approximations, this reduces to 
w* - w*(w,,) + w[2w w,, + (V2/c*)w .”] 
- (w.*w,,) = 0. (4) 


This is the basic equation of the theory presented 
here, but clearly approximate solutions are 
needed to clarify its meaning. 

First, if it be assumed that w.*<< w*, a con- 
dition easily met except for very low frequencies 
and quite near the earth, Eq. (4) simplifies to 


2. 2.0. 
w ww.) +2w w+ (V2 /c*)w 0 
The quantity 8w, /wz, occurring in the solution 


has magnitude 0.0043 and so can be ignored com- 
pared tol. Thus 


sa - 2 2)V/2 
we 2w [1 +(1 4Vw, /c*w yn. (5) 
This is precisely Gallet’s equation. 
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Next, if w<< wy, which is valid except far out 
and for quite high frequencies, Eq. (4) becomes 


w? - w[2w. + (V2w 2 /c*w,)|+ w” =0. 





Hence 
Vw 2 Vt 4 Vy fF) V2 
Dp p ps 
O90, *5y_ *\iew * ow . (6) 
H H H 


In the special case where w, << Vwy"/2c*wy, 
we find 
. 2 /p2 
weViw, /c*w 5» 

which is identical to Gallet’s equivalent result. 

Thus it has been shown that Gallet’s equations 
can be derived as limiting cases of a more 
general equation resulting from a quite different 
mechanism. 


(7) 


ee 


Several types of phenomena, such as the rising 
tones often associated with the “tails” of whis- 
tlers, are no doubt more easily understood in 
terms of the travelling-wave tube mechanism. 
However, it is likewise possible that the Doppler 
and proton-gyro concept can be made to cover 
cases not explainable by the former theory, and 
further work is progressing along these lines. 

Evidently, then, the two processes are com- 
plementary rather than mutually exclusive, and 
the full explanation of all types of dawn chorus 
will entail a combination of the two theories. 





* This research was generously supported by John 
B. Orr of Great Barrington, Massachusetts. 

'R. M. Gallet, Proc. Inst. Radio Engrs. 47, 211 
(1959). 

"See, for example, J. A. Ratcliffe, The Magneto- 
Ionic Theory and Its Applications to the Ionosphere 
(Cambridge University Press, Cambridge, 1959). 








SEARCH FOR AN ELECTRIC DIPOLE MOMENT OF THE ELECTRON* 


D. F. Nelson, A. A. Schupp, R. W. Pidd, and H. R. Crane 
The Harrison M. Randall Laboratory of Physics, The University of Michigan, Ann Arbor, Michigan 
(Received May 25, 1959) 


It has been pointed out'»? that the detection of 


an electric dipole moment (EDM) in an elementary 


particle would constitute a violation of both time- 
reversal invariance and space-reflection invar- 
iance. Experimental upper limits on the size of 
an EDM for several elementary particles have 
been found.*~’ A summary of these limits is 


given in Table I. Garwin and Lederman® have re- 


marked that the method used in searching for the 
muon EDM could be extended to higher accuracy 
for the electron. The anomalous magnetic dipole 
moment (MDM) experiment of Schupp, Pidd, and 


Crane® is an experiment of this general type. The 


Table I. Experimental upper limits of the electric 
dipole moments of elementary particles. 








Particle Electric dipole moment Reference 
Electron £107 om xe 3,4 
Positron £8107 om xe 4 
Proton £107 cm xe 5 
Neutron £3 x 107% cm xe 6 
Muon £2107" cm xe 7 
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purpose of this note is to show that this experi- 
ment does yield a new upper limit on the size of 
an EDM of the electron and to describe an exten- 
sion of the same experiment that is expected to 
have a considerably greater sensitivity for detect- 
ing an EDM of the electron. 

In the experiment on the anomalous MDM an 
unpolarized beam of electrons of 100-kev energy 
is sent parallel to a magnetic field. A portion of 
the beam that has been polarized by Mott scatter- 
ing from a thin gold foil at almost 90° is allowed 
to execute a helical path in the magnetic field (see 
Fig. 1). Application of a small trapping field 
allows the retention of the electrons in the nearly 
homogeneous magnetic field for times up to sev- 
eral hundred microseconds. The beam then 
strikes a second thin gold foil. Two Geiger 
counters, placed at a scattering angle of 80° and 
in orientations parallel and antiparallel to the 
magnetic field, are used to observe the Mott 
asymmetry as a function of the time the electrons 
are trapped in the field. The ratio of the counts 
in these two counters has a periodic variation 
with trapping time. The frequency of this varia- 
tion, the so-called difference frequency, is the 
magnitude of the angular velocity of spin pre- 








VoLUM 


POLARIZ 
TARGET 


FIG. 
the pol: 
indicat 
The pos 
Cy. 


cessiot! 
of an E 
quency 
part of 
the ele 
giving 
about t! 
frequer 
guishec 
anomal] 
quency 
gy (the 
ingly tc 
sensiti 
tained | 





- 





VYoLUME 2, NUMBER 12 PHYSICAL REVIEW LETTERS June 15, 1959 





—— 





The general equation” governing the precession 
of the spin vector 6 of an electron having an anom- 
alous MDM and an EDM in both a magnetic field 
B and an electric field E (which vary only ona 
macroscopic scale") is 

















do/dt =0 x6, (1) 
where 
a. ~~ = 5 Y VXE 
"‘ Moy E y+1 cc? 
2/5224 = iF 
Mo yY ov c 
gL [e- Bs aa. (2) 
FIG. 1. Experimental geometry. The precession of mo eae 
the polarization vector P about the magnetic field B is The factors a and f are defined in terms of the 
indicated at several positions along the beam trajectory. MDM vector 1 and the EDM vector p by p =(1+a) 
" positions of the two counters are shown by C, and xeh /2m,, and p =felic /2m,c; €, mo, and ¥ are 
the charge, rest mass, and velocity of the elec- 
tron; y=(1-87)"”* and B=v/c. The difference 
cession relative to the velocity. In the absence frequency is defined as 


of an EDM of the electron, this difference fre- 
quency is directly proportional to the anomalous ”D 
part of the MDM. The existence of an EDM of 
the electron would increase this frequency by 





whe OF (3) 


where We is the angular velocity of cyclotron 


> 


rotation. For the case of E =0 and v-B=0 


giving the electron spin an additional precession 





a 
about the vxB direction. The contribution to this Wy = w, (a? +£7B?)* , (4) 
frequency of the EDM can, however, be distin- 
guished by its energy dependence. Since in the where «, =eB/m,. (The error introduced by assum- 
anomalous MDM experiment the difference fre- ing v-B=0 is negligible.) It is seen that Wp in- 
quency was measured as a function of beam ener- creases with increasing energy and that at infi- 
gy (the magnetic field being changed correspond- nite energy wp =(a" + f?)”*, In the experiment on 
ingly to hold the beam at the same radius), a the anomalous MDM, wp*°*P >a when extrapolated 
sensitive test of the existence of an EDM is con- to infinite energy, but wy©®*P decreases with in- 
tained in the result. creasing energy. (See Fig. 2.) The most likely 
0.001166 
0.001165 
0.001164 5 
@ 0.001163+ 
3 
Sp ooonee } 
3 0.00116! 4 
re} 0.001160 + 
= i—t— THEORY, IN ABSENCE OF E.D.M. 
0.001159 + 
0.001158 + 
0.00157 + =a + + + + + + + 
0.0 02 0.4 0.6 0.8 1.0 L2 1.4 Lé Le 





vee (in 10° METERS® PER WEBER® ) 


FIG. 2. Plot of difference frequency data. The experimental values of the dif- 
ference frequency divided by eB/m, are extrapolated to infinite magnetic field 
(infinite energy). From the intercept (a*+f*)'”* the factor f (the EDM in units of 
eh /2m,c) can be found assuming the theoretical value of the g-factor anomaly a. 
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cause of this reversed slope is the presence of 

a small stray radial electric field. An electric 
field in any other direction, large enough to ex- 
plain the slope, would have to be about 100 times 
larger than one that could reasonably be expected 
to exist in the apparatus. If one represents the 
stray electric field as a small radial electric 
field E, constant in magnitude at the orbit radius, 
wp is given by 


2 2 4 
af] fea 
If a two-parameter (f and E) fit of Eq. (5) is made 


to the data, the theoretical value of a being as- 
sumed, the results are 


|p| =(1.6+0.6)x10~ cmxe, 
= - (0.048 +0.013) v/cm. (6) 


It is reasonable that a stray electric field of this 
size existed in the apparatus. The errors quoted 
are standard deviations of statistical error. 
Systematic errors are estimated to be of the or- 
der of 0.15% in each measurement of wp)/u». 
When these errors are included, the results are 
consistent with the nonexistence of an EDM of 
the electron. It is concluded that 


|p| <3x107"* cm xe. 


(7) 


It appears that the smallest EDM detectable by 
this method is about 3x10~"* cmxe. This would 
require more complete elimination of stray elec- 
tric fields (and measurement of them, if possible), 
improvement in the accuracy of the measurement 
of wp/w,, and extension of the measurements of 
wp/, to higher energy (~ 600 kev). 

The fundamental difficulty of this method of 
searching for an EDM is that the precession of 
the spin produced by the EDM is small compared 
with that produced by the anomalous part of the 
MDM. This difficulty can be removed if the mag- 
netic field in the system of the electron can be 
reduced enough so that the spin precession fre- 
quency due to the whole MDM equals the cyclo- 
tron rotation frequency. This can be done by 
providing a small radial electric field (~20 v/cm). 
In the absence of an EDM the effect of the radial 
electric field is to fix the spin vector relative to 
the velocity vector. If there is an EDM there will 
be a precession of the spin caused by the EDM 
around the radial direction (vxB), and there will 
not be any anomalous MDM precession present 
to obscure it. 

For the radial electric field scheme to work, 
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the polarization vector of the beam cannot initially 
be in the radial direction since there would then 
be no torque on the EDM. Since this is the initia) 
orientation of the polarization vector in the exper. 
iment on the anomalous MDM, a modification of 
this setup is necessary. The modification is to 
allow the polarization vector to precess, due to 
the anomalous MDM, to an orientation parallel 
to the velocity before switching on the radial 
electric field. The only further modification 
needed is the re-orientation of the counters in 
the radial direction. The ratio of the counts in 
the two counters will then vary with trapping 
time according to the frequency of the preces- 
sion produced by the EDM, given by 
Wy = Wf (8 +E/cB), (8) 

where the electric field is determined by 

a +(E/cB)[ (6 + ap? - 1)/8] =0. (9) 


The amplitude of this variation will be the same 
as that observed in the anomalous MDM experi- 
ment. A second pair of counters, oriented like 
those in the anomalous MDM experiment, will be 
used to determine experimentally the proper ad- 
justment of the radial electric field. The absence 
of a variation in the ratio of counts in these two 
counters when the trapping time is changed will 
indicate that the electric field is properly ad- 
justed. 

The practical limit of sensitivity of this method 
is imposed by the depolarization arising from the 
energy width of the beam in conjunction with the 
variation of the electric field with radius. Elec- 
trons of different energy will occupy orbits of 
different radii and so experience different elec- 
tric fields. At only one radius will the electric 
field be such as to eliminate the anomalous MDM 
precession exactly. At other radii the anomalous 
MDM will cause additional precession which will 
cause a depolarization. If a 1/r falloff of the 
electric field and an energy width of 100 ev at 
100-kev beam energy (caused by ionization losses 
in the first scattering foil) are assumed, a use- 
ful trapping time of 3 x10~* sec results. Since 
the phase of the EDM precession can be changed 
at will by changing the time the radial electric 
field is switched on, a tenth of a period of the 
EDM precession can be distinguished from spuri- 
ous effects. An upper limit on the EDM of the 
electron of about 4x10~'*® cm xe would result. 
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netic field, to stray electric fields, and to radial 
magnetic fields in the trapping region can be made 
negligibly small. 

We would like to thank Dr. L. M. Lederman for 
a stimulating communication. 





‘This work is supported by the U. S. Atomic Energy 
Commission. 
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‘These equations refer to expectation values of the 
various operators over a wave packet whose dimensions 
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vary appreciably. 





THEORY OF THE LATTICE VIBRATIONS OF GERMANIUM 


W. Cochran* 
Atomic Energy of Canada Limited, Chalk River, Ontario, Canada 
(Received May 15, 1959) 


The relation between frequency w and wave 
vector q (the dispersion relation) for the normal 
modes of germanium has been investigated ex- 
perimentally by Brockhouse and Iyengar’ (see 
also Brockhouse’) and by Ghose, Palevsky, 
Hughes, Pelah, and Eisenhauer.* The results 
do not agree with prior calculations made by 
Hsieh,* on the assumption that only nearest neigh- 
bor atoms interact and using the Born-von Kar- 
man theory of lattice dynamics. Nor is it suffi- 
cient to introduce second neighbor interaction’; 
indeed Herman’ has shown that force constants 
between atoms out to at least fifth neighbors are 
required to account for the experimental results. 
The theory then involves numerous parameters 
which have no clear physical significance. Lax® 
proposes to fit the data to a force model involving 
one parameter to represent electrostatic inter- 
action between quadrupoles generated by the lat- 
tice vibrations, and as many near neighbor param- 
eters as proves necessary. 

We have extended the Born-von Karman theory 
to apply to a model of the germanium crystal in 
which each atom is treated as a charged core 
coupled isotropically to an oppositely charged 
shell of negligible mass. This atom model has 
had some success in accounting for the dielectric 
properties’»® and lattice dynamics’ of alkali halide 
crystals. Its use introduces two long-range 
forces, firstly the electrostatic interaction be- 
tween dipoles formed by the relative displace- 
ment of cores and shells, and secondly a force 


which originates from the massless character of 
the shells. The shell model therefore simulates 
two effects which must occur in practice, electro- 
static interaction and relatively widespread re- 
distribution of valence electron density when the 
nuclei are displaced. The condition that the crys- 
tal as a whole can have no dipole moment® is auto- 
matically satisfied, so that the electrostatic in- 
teraction need not be treated as quadrupole- 
quadrupole. 

In the Born-von Karman theory the dispersion 
relation is determined by the condition for solu- 
bility of the equations 


pM, (kk’) - wm 


]U, (e’) =0, (1) 


poke’ 


where U(k) expi[q- F(Jk) - wt] is the displacement 
from fF(lk) of the kth atom in the /th cell. We de- 
fine My (kk’) by writing 


y= , ’ 
My (kk’) .., (kk oC. (kk’), 
where B,,(kk’) is related to the bonding force 


constant $,, (Ik, l’k’) between atoms (lk) and 
(l’k’) by the equation 


B , (kk’) =-)) 12, (Ue, Vk") 


x expiq: [F(1’k’) - F(7k)], (2) 


and C,,(kk’) is related to the electrostatic inter - 
action between charged atoms by a similar expres- 
sion.” For a diamond-type crystal with g along 
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[100] or [111] the set of Eqs. (1) can be reduced 
to three separate sets, each of which determines 
the dispersion relation for a single mode, and 
which can be written 


Di, AM(ek’) - wm, 8, ,,JU(e’)=0 (3) 


no kk’ 
where in general M is a linear combination of 
Myx; Myy, etc. For the shell model k =1...4, 
indices 1 and 2 referring to the two cores per 
unit cell, indices 3 and 4 to the corresponding 
shells, and we take m,=m,=m, m,=m,=0. The 





——__ 


to define new coefficients 
R = B(12) + B(32) + B(14) + B(34), 
T = B(32) + B(34), 


S =B(32). 
Each is a function of 4 and of two parameters. 


We introduce the abbreviations p=R/R,, 7 =T/T,, 


and o=S/S,, where (R)¢.,=-R,, and similarly 
for T, and S,. We also abbreviate 8 =a/v, where 
a@ is the atomic polarizability and v the volume 
of the (trigonal) unit cell. The dispersion rela- 
tion derived from the set of Eqs. (3) is then 








number of independent , (kk’) when only near- he el? A 
est neighbor units interact is six. It is convenient a =0, (4) 
A A, - mu? 
where 
A, =R, +B ~(1+B,+8C,)(1+ 171?) + 7(8C,* +B,o') + 1" (8C, +B,0) } (5) 
(1+B,+8C,)? - |BC,+B,o1l? 
* 
A=R+B -27(1+B, +8C,) +(8C, + Bo) + 7(8C," +Bo )t. () 
. (1+B,+8C,)? - |8C,+B,o1? 


Here B, and B, are constants, while C, and C, 
are dimensionless functions of q depending only 
on the geometry of the crystal structure. For 
example, 


, 9" 1 ate i 
Cy, - ar arm (EE), | , 


where f; is a lattice translation. These functions 
were evaluated by the Ewald method, as developed 
by Kellerman.” The fact that the elastic con- 
stants of germanium satisfy Born’s" identity now 
appears as a reflection of the fact that there are 
bonding interactions only between nearest neigh - 
bor units (cores or shells), and that p=7. Con- 
versely, the fact that the elastic constants of 
silicon and diamond do not satisfy the identity 
may merely indicate p#7, and not that there are 
other than first neighbor bonding interactions. 
Three of the seven parameters involved in the 
dispersion relation are determined by the elastic 
constants. On physical grounds we also set p =a, 
thereby eliminating a fourth parameter and mak- 
ing } y RR’) /®,(ek’) the same for kk’ =12, 14, 
and 34° The value of 8 deduced from the dielec- 
tric constant of germanium” is 0.100, leaving 
only two free parameters, the constants B, and 
B,. In fact, we also treated 8 as a disposable 
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parameter in fitting the data of Brockhouse and 
Iyengar, and obtained 8 =0.105, in satisfactory 
agreement with the value given by the dielectric 
constant. With the following numerical values, 


R, =2.923 x 105 dyne cm“, 
B, =4.099 x 10° dyne cm~, 
B, =3.200, 


B =0.105, 
and 


@ /& =0.6905 (in all three instances) 
xy xXx 


we obtained a good fit to the experimental data, 
as is shown in Table I and Fig. 1. The chief de- 








Table I, The elastic constants. 
Cy Cr, Cu 
Calculated 1.292 0.492 0.676 
Measured, 300°K 1. 288 0.483 0.671 
Measured,> 0°K 1.310 0.490 0.687 








*H. J. McSkimin, J. Appl. Phys. 24, 988 (1953). 
bm. E. Fine, J. Appl. Phys. 26, 862 (1955). 
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FIG. 1. A comparison of the theoretical dispersion 
curves with the experimental results of Brockhouse 
and Iyengar. 


fect is that the calculated frequency of the LA 
mode is too low when q is near the zone boundary 
inthe [111] direction. 

It is remarkable that such a simple model should 
fit the data as well as it does. Our procedure, 
however, finds theoretical backing in the work of 
Mashkevich and Tolpygo.** Using the Hartree- 
Fock and adiabatic approximations, they have 
shown that the energy perturbation should be ex- 
pressed as a quadratic function of the nuclear dis- 
placements and the atomic dipole moments. It 


can be shown that the shell model represents the 
most general (and probably the simplest) way of 
accomplishing this. 

I am greatly indebted to Dr. B. N. Brockhouse 
for many helpful discussions and suggestions. 





* 

On leave from Crystallographic Laboratory, Caven- 
dish Laboratory, Cambridge, England. 
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DETECTION OF COLLOIDAL CENTERS IN LITHIUM HYDRIDE BY ELECTRON RESONANCE 


W. T. Doyle,* D. J. E. Ingram, and M. J. A. Smith 
Department of Electronics, Southampton University, Southampton, England 
(Received May 18, 1959) 


This Letter reports electron resonance meas- 
wements on lithium hydride which has been ir- 
radiated with ultraviolet light and the results 
show that the lithium hydride contains colloidal 
tenters of lithium metal, with diameters small 
tompared with the microwave skin depth. The 
lithium hydride was irradiated for about fifteen 
ninutes with 2537A radiation from a 250-watt 
mercury lamp and a fairly dark coloration was 
produced. These samples were then inserted in 
13-cm wavelength electron resonance spectrom- 
tter employing 100-kc/sec field modulation and 
thase- sensitive detection, and their electron 


resonance spectra recorded. All the irradiated 
samples gave a single resonance line at g=2.002 
+0.001 with a half-width of 0.3 gauss—the latter 
being determined by the magnetic field inhomo- 
geneity. In contrast to this, no signals at all 
were obtained from the unirradiated material. 
This very narrow signal shows that the para- 
magnetic entity cannot be due to normal electron 
excess centers, which have wide inhomogeneously 
broadened lines.’ It was therefore postulated 
that the absorption was in fact due to the conduc- 
tion electrons in colloidal particles of lithium. 
In order to test this hypothesis a careful inves- 
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FIG. 1. Electron spin resonance of irradiated LiH. 
Curves (a) and (b) show the derivative of the absorption 
signals obtained from LiH and from a standard free- 
radical sample, respectively, at room temperature. 
Curves (c) and (d) show the corresponding signals at 
liquid oxygen temperature. 


tigation on the temperature dependence of the 
intensity of the electron resonance was made, 
and the results are summarized in Fig. 1. 
Measurements were made at room temperature 
and at 90°K, and in each case the curves from 
the lithium sample and a standard free-radical 
sample were recorded under identical conditions. 
It can be seen from the figure that, whereas the 
intensity of the signal from the free radical in- 
creases markedly as the temperature is lowered, 
as would be expected, that from the irradiated 
lithium hydride remains constant. The fact that 
narrow symmetrical electron resonance lines 
are obtained with intensities independent of tem- 
perature shows that they must originate from 
colloidal particles of lithium with diameters 
small compared with the skin depth at micro- 
wave frequencies.’ 

Further confirmation of the presence of such 
colloidal particles is obtained from the optical 
absorption spectra. Optical absorption meas- 
urements have been made on a very small single 
crystal and the differential absorption of the ir- 
radiated and nonirradiated crystal is plotted in 
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FIG. 2. Differential optical absorption spectrum of 
irradiated LiH. Optical density is plotted against photon 
energy in ev. 


Fig. 2. It is evident that the irradiation has pro- 
duced a strong band at 650 my. The position of 
this band does not agree at all with the positions 
predicted for F or aggregate bands in the hy- 
drides,* whereas a strong band at 2512 A has 
been reported‘ and its position agrees well with 
that predicted for the F band. Preliminary esti- 
mates on the expected position of the colloid ab- 
sorption band, based on the index of refraction 
of the host medium® and the plasma wavelength® 
of the alkali metal, can be made, however. 
These preliminary calculations suggest that, 
although very small colloidal particles would be 
expected to give rise to an absorption at a con- 
siderably shorter wavelength than 650 my, it is 
nevertheless quite possible for the band to be 
shifted to this region if larger colloidal particles 
are involved. The electron resonance measure- 
ments show that at least 10** unpaired electrons 
per gram are present in the specimens, and 
hence a relatively large colloidal content of lith- 
ium is present. 

The general conclusion is that both the electron 
resonance and optical absorption spectra show 
that colloidal centers of metallic lithium are 
formed when lithium hydride is irradiated with 
ultraviolet light. It would appear that this may 
be a very general effect, and initial studies on 
the thermal and optical decomposition of various 
azides have also yielded narrow electron spin 
resonance lines. Further work on these and 
other centers, employing much more homogene- 
ous magnetic fields, is now in hand so that the 
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theories relating the width of the lines, the spin- 
orbit coupling, and the mean free path of the con- 
duction electrons, can be checked in detail. 
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National Science Foundation postdoctoral fellow. On 
leave from Dartmouth College, Hanover, New Hamp- 
shire. 
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EFFECTIVE EXCHANGE CONSTANT IN YTTRIUM IRON GARNET* 


D. T. Edmonds and R. G. Petersen 
Department of Physics, University of California, Berkeley, California 
(Received May 18, 1959) 


We have measured the specific heat of a pure 
(99.99 % yttrium) polycrystalline sample of yt- 
trium iron garnet (YIG) in the temperature range 
of 1°K to 4°K. We find that the specific heat may 
be analyzed into a part that depends on tempera- 
ture as T*” (the spin-wave contribution) and a 
part that depends on temperature as T® (the 
contribution from the lattice) as illustrated in 
Fig. 1. 

We thus obtain 


C = 68.97% +29.3T° erg/cm® deg K, (1) 


where C, is the specific heat per cm*® based on 





T T T T T T T T q 











0 ! | | | | ! ! | | 
4 5 6 8 9 10 
Te [x] 

FIG. 1. The specific heat per unit volume of poly- 
crystalline YIG analyzed into its two components by 
plotting C,,/T*/* against T*’?. The points give the re- 
sults obtained in one experiment. Another experiment 
performed on the same specimen at the extremes of 
the temperature range yielded results in agreement 
with those plotted to within 2%. 





the x-ray density’ of 5.19 g/cm’. 
Using the equation 


C (lattice) = (12/5) n* Nk ,(T/8)°, (2) 


where kp is the Boltzmann constant, we obtain 
for the Debye temperature 6 =454°K, if we take 
N to be the total number of atoms of all kinds in 
unit volume. 

Assuming a dispersion relation for the spin 
waves of 


fiw =Dk*, k=2n/x, (3) 
one has the standard result for the magnon spe- 
cific heat per unit volume: 
C,, (magnons) =[(1 5/32)(1.341)0-** ke (ek ,7/D)™ 
- 2 
=0.113 k p(k ,T/D) : (4) 


The constant D may be expressed in terms of 
the Landau-Lifshitz exchange constant A, the 
saturation magnetization M,, and the spectro- 
scopic splitting factor g as 


D=2Ag(u,/M,), (5) 
where pp is the Bohr magneton. In this manner 


we obtain from our measurement of the spin- 
wave contribution to the specific heat, 


A(YIG) = 0.192 x10~* erg/cm, (6) 


where we have taken the value® 41M, =1740 gauss 
for YIG. This is to be compared with a value of 


A(Fe) =2x10°* erg/cm 


in iron. If we analyze in this manner the results 
of Kouvel* who measured the specific heat of 
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magnetite in the same temperature range, we 
obtain 


A(Fe,O,) = 0.232 x10-* erg/cm. 


Suhl® has shown that in a sphere, z-directed 
spin waves having a wave number given by k, in 


he = 20M .*/3A 


are degenerate in frequency with the uniform 
precessional mode (k=0). We expect*»” maxi- 
mum coupling of the uniform precessional mode 
to transverse phonons via the z -directed spin 
waves when the frequency is given by 


i =v,k,/2n, 


where v; is the velocity of transverse phonons. 
Using our value of A and the value v,=3.87 «10° 
cm/sec obtained in polycrystalline YIG at 25°C 
by McSkimin,® we find 


|Ry| = 4.60 x10° cm™ 


and 


f= 2.83 x10'° sec™* or d= 1.05 cm. 


It is our pleasure here to record our debt to 
Professor C. Kittel who suggested the investiga- 
tion and to M. K. Jack of Hughes Research Labora- 
tories who kindly supplied us with the specimen 
of YIG. 





_ 

Supported in part by the National Science Foundation, 

'F. Bertaut and F. Forrat, Compt. rend. 242, 382 
(1956). 

*C. Herring and C. Kittel, Phys. Rev. 81, 869 (1951), 

3w. P. Wolf and G. P. Rodrigue, J. Appl. Phys. 29, 
105 (1958). > 

‘J. S. Kouvel, Phys. Rev. 102, 1489 (1956). 

5H. Suhl, Proc. Inst. Radio Engrs. 44, 1270 (1956); 
J. Chem. Phys. Solids 1, 1 (1957). 

SC. Kittel, Phys. Rev. 110, 836 (1958). 

"Akhiezer, Bar’iakhtar, and Peletminskii, J. Exptl. 
Theoret. Phys. U.S.S.R. 35, 228 (1958)[translation: 
Soviet Phys. JETP 8, 157 (1959)]. 

8H. J. McSkimin (private communication). 





SUPERCONDUCTIVITY VERSUS FERROMAGNETISM IN LANTHANUM-GADOLINIUM ALLOYS 


R. A. Hein and R. L. Falge, Jr. 
United States Naval Research Laboratory, Washington, D. C. 


B. T. Matthias and C. Corenzwit 
Bell Telephone Laboratories, Murray Hill, New Jersey 
(Received May 25, 1959) 


The differential paramagnetic susceptibility of 
samples of LaGd alloys containing from 0.60 to 
5.00 % Gd has been measured in the temperature 
range from 4.22 to 0.10°K. Matthias et al.’ have 
shown that the addition of small amounts of Gd to 
La produces solid solutions with unusual mag- 
netic properties. Alloys containing up to 0.90% 
Gd are superconductors with transition tempera- 
tures (Te) which decrease linearly with increas- 
ing Gd content. Alloys containing 3.0 or more 
percent Gd are ferromagnetic with Curie tem- 
peratures which increase with increasing Gd con- 
tent. The alloy containing 3.0% Gd becomes fer - 
romagnetic at 1.3°K, the lowest temperature they 
could attain. Extrapolation of these data suggests 
that the T,. vs percent Gd curve and the Curie 
temperature vs percent Gd curve would intersect 
for an alloy of 1.25% at a temperature of 0.5°K. 
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Thus two possibilities are suggested: (1) the 
curves meet and terminate, indicating the pos- 
sibility of a ferromagnetic superconductor for 
one specific composition, or (2) they intersect 
in such a way that there would exist a range of 
compositions for which a given alloy would ex- 
hibit (in zero magnetic field) both superconduc- 
tivity and ferromagnetism, but at different tem- 
peratures. The primary purpose of the work re- 
ported here was to investigate the validity of 
these extrapolations. 

Temperatures below 1°K were produced by the 
magnetic cooling method. The sample to be in- 
vestigated was cemented in a copper holder which 
was in thermal contact with the paramagnetic 
salt via a copper rod 12 cm long and 3 mm in 
diameter. A carbon resistor cemented directly 
on the sample indicated that it was in intimate 
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FIG. 1. Galvanometer deflection as a function of 
temperature obtained in zero applied magnetic field for 
the alloys containing 1.25 and 1.00% Gd. 


thermal contact with the salt.? The differential 
magnetic susceptibilities of the salt and sample 
were independently observed by means of a dc 
mutual inductance technique described in detail 
elsewhere.® 

In Fig. 1 galvanometer deflections obtained in 
ero applied magnetic field with the 1.25 and 
1.00 % Gd alloys are plotted as a function of tem- 
perature. Such a plot represents susceptibility 
inarbitrary units vs temperature. These curves 
suggest the occurrence of a magnetic transition 
of some kind. 

One possible interpretation of the curves in 
Fig. 1 is that the decrease in the susceptibility 
at the lower temperatures is due to supercon- 
ductivity in parts of the sample. The entire sam- 
ple could not be superconducting as the suscepti- 
bility was still paramagnetic at the lowest tem- 
peratures. If the decrease in susceptibility were 
caused by the diamagnetism of superconducting 
inclusions, the application of a magnetic field 
less than critical would have no effect on the sus- 
teptibility while a field greater than critical 
would cause an increase in the susceptibility. 
Application of a magnetic field actually caused 


the susceptibility to decrease, as shown in Fig. 2. 


Therefore, these alloys (1.00, 1.10, and 1.25% 
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FIG. 2. Galvanometer deflection as a function of 
temperature and magnetic field for the alloy contain- 
ing 1.10% Gd. 


Gd) do not become superconducting in the earth’s 
magnetic field down to a temperature of 0.1°K.* 

The maxima in the curves of Figs. 1 and 2 can 
be interpreted as evidence of antiferromagnetic 
or ferromagnetic transitions. Although pure Gd 
is ferromagnetic, Thoburn et al.,° on the basis 
of magnetic moment measurements at low tem- 
peratures, reported antiferromagnetic transitions 
in LaGd alloys containing 47 to 90% Gd. On the 
other hand, Matthias et al.’ observed large rem- 
anent moments for LaGd a alloys containing 3 to 
10 % Gd, and concluded that these alloys are fer- 
romagnetic below 6.0°K.’ In keeping with this 
latter conclusion the temperatures of the maxima 
in the susceptibility have been taken, in this re- 
port, to be approximations to ferromagnetic 
Curie temperatures.* Curie temperatures so de- 
fined are plotted as a function of Gd content in 
Fig. 3. 

The 0.60 % Gd alloy becomes superconducting 
at 2.6°K.' Extrapolation of the Curie tempera- 
ture curve in Fig. 3 indicates that this alloy 
might also become ferromagnetic at 0.3°K. The 
susceptibility of a sample of this alloy was para- 
magnetic at 4.2°K and increased as the tempera- 
ture was lowered. A superconducting transition 
was observed which, in zero applied magnetic 
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FIG. 3. Superconducting transition temperatures and 
ferromagnetic Curie temperatures of the La-Gd alloys 
as a function of Gd concentration. 


field, was spread over the temperature interval 
from 2.8°K to 1.3°K. From 1.3°K down to 0.1°K 
the sample remained diamagnetic with a constant 
susceptibility. No evidence of ferromagnetism 
was observed in zero magnetic field. 

Alloys in the range of 0.80 to 0.98% Gd were 
also studied. All samples were found to be super- 
conducting. The zero field transitions (7,) for 
the 0.60 and 0.80 % Gd samples were in accord 
with the linear behavior previously reported.’ 
However, departure from linear behavior was 
noted for a 0.95% sample. Its T, was higher 
than expected. The first reaction to this result 
was to question the Gd content of the sample. 
Further doubt about composition arose when a 
0.98% sample was found to possess a 7,, which 
agreed with the linear dependence. To resolve 
this question of composition, a second 0.95% 
sample was fabricated. Its T,, was higher than 
that of the first. In fact, its T,, was even higher 
than that of the 0.80% sample. A 0.90% sample 
was then investigated and the results substanti- 
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ated the deviation from linearity observed for the 
0.95% samples. The dependence of T,, on Gd 
concentration is depicted in Fig. 3. In drawing 
the curve in this figure, more weight was given 
to the second 0.95% sample as extra care was 
taken in its fabrication. 

In Fig. 3 one notes rather strange behavior in 
the vicinity of 0.95% Gd. Increasing the Gd con- 
tent from 0.95 to 0.98% causes T,, to decrease 
by a factor of two. Increasing the Gd content 
from 0.98 to 1.00 % causes the sample to exhibit 
ferromagnetism rather than superconductivity. 
On the other hand, decreasing the Gd content 
from 0.95 to 0.80 % has relatively little effect on 
T,. The sharp decrease in T,, observed for those 
alloys between 0.95 and 0.98 % is consistent with 
the ideas of Suhl and Matthias.’ They predicted 
that the slope of T,, vs percent Gd would approach 
infinity as T, approaches zero. The observed 
slope is 30°K per percent Gd. The behavior of 
T¢ in the region between 0.80 and 0.95 % and the 
change from superconductivity to ferromagnetism 
between 0.98 and 1.00 {are not yet understood. 

Measurements of the susceptibility of the 0.954 
sample in magnetic fields sufficient to destroy 
superconductivity suggest that this sample is fer- 
romagnetic at the fields and temperatures involv 
Critical magnetic field measurements, as well as 
magnetic moment measurements, will be includ 
in a more comprehensive report. 
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LONGITUDINAL MAGNETORESISTANCE OF METALS IN HIGH FIELDS. 


B. Lithi 
Institut fur Kalorische Apparate und Kaltetechnik, Swiss Federal Institute of Technology, Zurich, Switzerland 
(Received May 4, 1959) 


Recently de Launay, Dolecek, and Webber’ have 
measured the longitudinal resistance change of 
polycrystalline copper in magnetic fields up to 
100 kilogauss, and they found a saturation al- 
though the transverse magnetoresistance appears 
to increase linearly without limit. We want to 
present, here, preliminary results of longitudi- 
nal magnetoresistance measurements in high 
magnetic fields on a number of other polycrys- 
talline metals which show a similar effect. 

Magnetic fields up to 200 kilogauss were pro- 
duced by condenser discharges through coils 
cooled by liquid air or liquid helium.” Pulse 
lengths between 5 msec and 30 msec were avail- 
able.* A number of experimental problems must 
be solved for such measurements. Shaking of the 
specimen due to shaking of the field coil and due 
to the Lorentz forces was prevented by mounting 
the specimen very carefully; the eddy-current 
heating was diminished by taking wires with only 


“3 0.2-mm diameter, and by choosing an appropriate 


discharge time. Furthermore, care had to be 


487 taken to place the specimen parallel to the exter- 


nal field, since the transverse magnetoresistance 
can be orders of magnitude bigger than the longi- 
tudinal one. The ratios of the transverse to the 
longitudinal magnetoresistance in the fields attai- 
table with our apparatus are, for example, for Li, 
Ag, Au, and Pt about 10, for Sn 10*, and for Pb 
and Zn about 10°. As a control of our technique, 

the results on some metals were checked by 

means of dc measurements in low fields on speci- 
mens of extreme purity. 

Figure 1 shows the results of Li, Ag, Au, Zn, 
in, Ni, Pt, and Pb in a Justi-Kohler plot, in 
vhich AR means the resistance change, R,, the 
ero-field resistance at temperature T (liquid air 
or liquid helium) and R, the resistance at the 
Debye temperature. The results plotted in this 
figure show quite clearly that there exists a very 
strong trend to saturation in the high-field region 
for all metals measured so far except Li, which 
shows a strict linear behavior, and Pb, which 
llows an H“? law. Bearing in mind the above- 
nentioned experimental difficulties, one may 
state that the longitudinal magnetoresistance of a 
great number of polycrystalline metals shows a 
very weak field dependence in high magnetic fields. 


pi 
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FIG. 1. Longitudinal magnetoresistance of 
@ Li (430) caAg (220) @Au (320) 
A Zn (240) wSn (160) O©ONi (410) 
V Pt (240) 4 Pb (90). 


The figures in brackets indicate the Debye 6 values 
used in calculating Rp. 


For single crystals this behavior would be in 
nice agreement with existing theories.‘ How- 
ever, as Ziman’ has pointed out, in the case of 
polycrystalline metals one has to take an average 
over all microcrystal directions, and here the 
special form of the energy surface plays an im- 
portant role. Ziman’s model with some very 
special assumptions gives a linear dependence 
for the transverse and the longitudinal magneto- 
resistance, but this model with its cylindrical 
surface is too special to account for a real metal 
such as Li. For the other monovalent metals 
Cu, Ag, and Au, one can give an explanation of 
the magnetoresistance behavior (linear depend- 
ence for the transverse case, saturation for the 
longitudinal one) by allowing the Fermi surface 
to touch the Brillouin-zone face. However, 
since our results show that a saturation-like be- 
havior of the longitudinal magnetoresistance 
appears to occur for a wide range of metals and 
with only few exceptions, an explanation inde- 
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pendent of the detailed form of the energy sur- 
faces is needed. These measurements are being 
continued on other substances. 

I wish to thank Professor P. Grassmann, the 
Director of this Institute, and Dr. J. L. Olsen 
for their kind encouragement and support. This 
work was supported financially by a Arbeitsbes- 
chaffungskredit des Bundes. 
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STIMULATED EMISSION OF RADIATION BY RELATIVISTIC ELECTRONS IN A MAGNETIC FIELD’ 


Jurgen Schneider 
Department of Physics, Duke University, Durham, North Carolina 
(Received April 27, 1959) 


The Schrédinger equation of an electron, mov- 
ing perpendicularly to a magnetic field H, can be 


reduced to the equation of the harmonic oscillator.’ 


The kinetic energy W; of an electron in the state 
i is therefore given by W; =(i+3)iuw,, where 

w, =eH/m,c is the cyclotron frequency, associ- 
ated with all transitions between the states i and 
i+1. Using the relativistic Schrodinger equation 
and neglecting the spin, the kinetic energy levels 
of the electron are found to be? 


1 
W, =m,c*[1 + 2(i +4)iw, /m,c*}° -m,c*. (1) 
Transitions between states i+1 and i now result 
in a whole spectrum of frequencies w; , 1, ;: 


ars @,(1 - iw, /m,c’). 


(2) 


Second order terms in fw, /m,c? have been neg- 
lected. Consider N; electrons in the energy state 
i. Under the influence of an alternating electric 
field they will undergo absorption and induced 
emission at the frequencies w, 1 and w; 
The net transfer of power is 


P; =N Pe; 


i,i+ 4¢<2° 


ii +1" %,i+1- nail ©) 
The transition probability w; ;,1=wj41,;1 


given by 


=A Ey i41) 8). (4) 


”3,i+1 
The quantity 


a 
= ; 2 
Meet =e[(i+1)R/2m, a, | 


(5) 


is the matrix element of the harmonic oscillator 
for the transition between the state i and i+1; 
E is the amplitude and w the frequency of the ex- 
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ternal electric field. For sufficiently sharp lines, 
i.e., for a small collision frequency 1/7, the re- 
sponse curve g(w) can be assumed to be Lorent- 
zian: 


Bw) =T/A1+(w, sy - w)?7?]. (6) 


Introducing Eqs. (2), (4), (5), and (6) into Eq. (3), 
one obtains 
P,=N (e*/my)(E*/2)¢, 


(1 - ia) (i +1) [1-G-a} =) (" 


= = 2 - 1 
1+(w, 544 w)*7? 1+(w, a w)? 





where a =hw,/m,c’. Neglecting again second 
order terms in a, the function ¢ can be written 


W x 
g(x) = at - +x? ee (1 +x")? , (8) 








= > aS 
-0.5 \/ ° +0.5 
xX —e 
FIG. 1. Cyclotron resonance absorption of a non- 


relativistic (A) and a relativistic (B) electron. 
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where x =(w; ;,1~-)7, W is the kinetic energy 
of the electron, and Q=w,7. In the nonrelativistic 
limit, i.e., for W/m,c* =0, the cyclotron reson- 
ance line is Lorentzian. However, for QW/m,c*>1, 
the function g may become negative, corresponding 
to a net stimulated emission instead of absorption. 
This effect is shown in Fig. 1 where we have taken 
QW/m,c? =1.5, e.g., Q=10000 and W =76 ev. 

If the electrons are not monochromatic, but 
have rather a distribution of energies, one ob- 
tains essentially the same effect, if there is an 
overpopulation of the upper states.” The above 
formulas (7) and (8) have also been derived? by 
means of the Boltzmann transfer equation, tak- 


ing into account the dependence of the electron 
mass on the kinetic energy. 

It does not appear unlikely that this effect could 
be used for a new type of maser, which would re- 
quire no microwave “pump” and no low-tempera- 
ture operation.’ 

The author wishes to thank Dr. Walter Gordy 
for his interest in this work. 





7 

This research was supported by the U. S. Air Force 
through the Air Force Office of Scientific Research of 
the Air Research and Development Command. 

'L. Landau, Z. Physik 64, 629 (1930). 

23. Schneider (to be published). 





STRENGTHS OF HIGH-ENERGY CAPTURE GAMMA RAYS IN W'* 


D. J. Hughes, M. K. Brussel, J. D. Fox, and R. L. Zimmerman 
Brookhaven National Laboratory, Upton, New York 
(Received May 27, 1959) 


During the past few years the Brookhaven fast 
chopper’ has been devoted to an extensive study 
of the parameters of energy levels in heavy 
nuclides excited by the capture of slow neutrons. 
The results, obtained primarily from total cross 
sections alone, have given information on the dis- 
tribution laws of radiation widths, neutron widths, 
and spacings of these energy levels.” However, 
little is known of the spins, J, of the levels, 
aside from the limitation to two values, /+ 4, 
with J the target nucleus spin. In a few favorable 
cases, for which I, is larger than ry, the total 
cross section suffices for J determination, but 
in general, partial cross sections, scattering or 
capture, must be measured as well. 

Another possible method for spin measurement 
is based on the capture gamma rays emitted from 
the individual levels. The most direct use of cap- 
ture gamma rays for determination of J applies 
when the ground-state transition is allowed for 
one of the spins only.* An example‘ is the com- 
pound nucleus W’**, in which the electric dipole 
ground-state transition, of energy’ 7.42 Mev, is 
allowed for the J=1 but forbidden for the J=0 
levels. In this ground-state technique of spin 
identification the question of the variation in 
transition probability from level to level is im- 
portant, for it might be possible that a J=1 level 
would exhibit no observable ground-state gamma 
because of an extremely weak transition proba- 
bility. Because of this application to spin deter - 


mination, and the intrinsic interest in transition 
probabilities as well, it is desirable to gain in- 
formation on the high-energy capture gamma 
rays in heavy nuclides. 

Essentially no information exists concerning 
the relative transition probabilities involved in 
the present situation, that is, the emission of 
gamma rays from neighboring levels a few ev 
apart to the same final state, some 6 or 7 Mev 
distant. The typical neutron capture gamma-ray 
measurements do not supply the necessary in- 
formation, for they are limited to excitation by 
one energy only, usually thermal neutrons. Elec- 
tric dipole transitions in the hundred-kev range® 
show wide variations in transition probability and 
are usually hindered by large amounts relative to 
the high-energy capture gamma rays; the latter 
are in reasonable agreement with estimates® based 
on Weisskopf’s single-particle model, although 
exhibiting variations of the order of a factor of 
ten from nuclide to nuclide. 

Earlier capture gamma measurements with the 
Brookhaven fast chopper showed‘ strong ground- 
state transitions’ for several of the known res- 
onances in the target nucleus W'®*, indicating that 
J=1 for these resonances. In that work gamma- 
ray energies were measured in sodium iodide 
crystals located near the capturing sample, and 
neutron energies simultaneously by time of flight 
to the sample. The same equipment has now been 
installed at the NRU reactor, Chalk River, Can- 
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ada, where the fast chopper flux is some 80 times 
greater and correspondingly higher resolution is 
possible. Definite ground-state gamma rays have 
now been detected in five of the six known levels 
up to an energy of 100 ev; thus these levels are 
identified as J=1. The spins are in agreement 
with the results of Waters et al.* based on par- 
tial cross sections, except for the 41- and 66-ev 
levels, too weak to be analyzed by partial cross- 
section methods. It is noteworthy that ground- 
state gammas can be detected by the present 
method for such weak levels. 

The transition probabilities of the ground-state 
gamma rays for the five J=1 levels are given in 
Table I relative to the average value, taken as 
unity. The absolute value of this average proba- 
bility, obtained by comparison with thermal 
neutron results, *»*® is 2.6x10~* ev, expressed as 
a partial radiation width. The striking fact about 
the transition probabilities is their small varia- 
tion from level to level, an average deviation 
from the mean of only 20% which is very different 
from the wide variation of neutron widths, also 
listed in Table I for comparison purposes. These 
results are surprising in the light of the discus- 
sion of Porter and Thomas,” who assume that the 
individual gamma rays would show an extreme 
variation from level to level of the same type as 
that observed" for neutron widths. Although the 
number of levels in the present work is small 


Table I. Relative transition probabilities of ground- 
state gamma rays emitted from specific J=1 levels 
in W'*4, The levels are identified by the energy of the 
corresponding neutron resonance and the reduced neu- 
tron widths of these resonances are also listed. 








Reduced Relative 
Neutron neutron transition 
energy widths probability 
(ev) (1075 ev) (mean value = 1. 00) 
7.6 1.2 1.27+0.24 
27.1 9.2 1.01 20.03 
41 0.16 0.53 + 0.12 
47 37 1.00 + 0.04 
66 0.26 1.19 + 0.27 
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as yet, the close grouping of the observed transi- 
tion probabilities is suggestive that the high- 
energy radiation from neighboring levels in heavy 
nuclides does not exhibit a range as wide as the 
Porter -Thomas neutron width distribution. Meas- 
urements for other similar nuclides are continu- 
ing because of the intrinsic interest of the transi- 
tion probabilities and for the determination of 
spins of more levels.’” We wish to acknowledge 
the cooperation of the Chalk River scientists in 
setting up the cooperative program that made 
these results possible, particularly H. B. Lewis 
and C. H. Westcott. H. Palevsky and R. E. Chrien 
have given valuable help in the design of the equip- 
ment and W. C. Olsen in its operation. 





‘Seidl, Hughes, Palevsky, Levin, Kato, and Sjéstrani, 
Phys. Rev. 95, 476 (1954). 

?For a summary see D. J. Hughes, Neutron Cross 
Sections (Pergamon Press, New York, 1957). 

H. H. Landon and R. Rae, Phys. Rev. 107, 1333 
(1957). 

‘Fox, Zimmerman, Hughes, Palevsky, Brussel, and 
Chrien, Phys. Rev. 110, 1472 (1958). 

5B. B. Kinsey and G. A. Bartholomew, Can. J. Phys, 
31, 1051 (1953). 

*B. B. Kinsey, Handbuch der Physik (Springer-Verlag, 
Berlin, 1957), Vol. 40. 

"These might include the similar allowed transitions 
to the 2+ state at 111 kev, but none of the present con- 
clusions would be significantly affected. 

‘Waters, Evans, Kinsey, and Williams, Nuclear Phys. 
(to be published). 

*Groshev, Demidov, Lutsenko, and Plevhanov, Atlas 
of Capture Gamma Rays (Moscow, 1958). 

0C. E. Porter and R. G. Thomas, Phys. Rev. 104, 
483 (1956). 

‘Dp, J. Hughes and J. A. Harvey, Phys. Rev. 99, 
1032 (1955). 

12Similar measurements of ground-state transitions 
for the same five levels in W'*‘ have recently been re- 
ported by Huynk, Julien, Corge, Netter, and Simic, 
Compt. rend. 248, 2330 (1959). Their results and 
those of Table I are in excellent agreement. In addi- 
tion, J. R. Bird has reported qualitative results for 
these levels in the Proceedings of the Second United 
Nations International Conference on the Peaceful Uses 
of Atomic Energy, Geneva, 1958 (United Nations, Gen- 
eva, 1959), paper P/35. 
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SYSTEM OF LEVELS IN EVEN-EVEN NUCLEI" 


C. A. Mallmann 
Argonne National Laboratory, Lemont, Illinois 
(Received May 11, 1959) 


It is of interest, in order to stimulate further 
experimental and theoretical work, to point out 
the existence of a system of levels in even-even 
nuclei whose energies are related by an empiri- 
cally determined law which seems to be appli- 
cable throughout the nuclear chart. 

If one plots the ratio R of the energy of the 
first excited 6+ (denoted '6+) and 8+ ('8+) states 
to the energy of the first excited 2+ ('2+) state as 
a function of the ratio R('4+) of the energy of the 
first excited 4+ ('4+) state to the energy of the 
first excited 2+ ('2+) state, one finds that prac- 
tically all the points lie on two curves, one for 
the ‘6+ states and one (tentative) for the '8+ 
states; see Fig. 1. Some of the energy ratio 
[for R(*4+) larger than 2.9] were taken from 
Scheuer and Aisenberg’s paper’ with the excep- 
tion of Gd*®*.* All the values used are given in 
Table I. Question marks in Fig. 1 mean that 
there is uncertainty in the assignment of the levels 
considered. We also show in Fig. 1 the obvious 
lines for R('4+) and R('2+) =1 in order to visualize 
better the position of the '2+, 14+, *6+, *8+ levels. 
Naturally the '0+ ground state is at R =0. 

As one would expect, the R('8+) and R('6+) 
curves coincide over a range of R('4+) values with 
the straight lines R('8+) = -312/7 + (594/35) R(#4+) 
and R(?6+) = -11+(27/5)R(?4+). These are a di- 
rect consequence of 


E = al(I+1) - b+), (1) 
which is the expression for the energy of a rota- 
tional state in an axially symmetric nucleus* tak- 
ing into account the rotational -vibrational‘* inter - 
action, where a and b are constants and 7 the spin 
of the state. The coefficients in the expressions 
for R(*4+) and R(*6+) are independent of a and b. 
The curves show that the rotational nuclei are 
confined to R(*4+) values larger than or equal to 
3.27 because for values smaller than this the 
empirical curves deviate from (1). Of the rota- 
tional nuclei, Hf'”*, U*°*, U?%*, and Pu*® do not 
fall on the theoretical curve (1) within the experi- 
mental errors. 

The theory of Davydov and Filippov’ which treats 
rotations of nonaxially-symmetric nuclei, extends 
the agreement between theoretical and empirical 


curves down to values of R('4+) =8/3 because 
R('4+) and R('6+) are reduced from 10/3 and 7, 
respectively, for y =0° to 8/3 and ~5.0 for y =30°. 
Unfortunately the precise predictions for R('6+) 
have not been published so that a more careful 
comparison with experimental values cannot be 
made. The R(?2+) (second 2+ excited state) values 
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FIG. 1. The figure shows the ratios R(!g+), R('6+), 
R('4+), and R('2+) as a function of the ratio R('4+). 
Experimental points are indicated with a circle for 
R('8+) and a plus sign for R('6+). The R('6+) points 
for R('4+) values in the range from 3.28 to 3.33 are 
shown in a magnified scale at the upper-left side of 
the figure. Experimental errors are not indicated in 
order to avoid confusion in the graph. They are given 
in Table I. Question marks mean that there is an un- 
certainty in the assignment of the levels considered. 
Full lines are experimental; the dashed line is tenta- 
tive for R('8+). The dotted lines are the theoretical 
predictions for rotational states in axially symmetric 
nuclei taking into account the rotational-vibrational 
interaction. 
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Table I. Experimental energy ratios with their errors used in Fig. 1. 
E('2+ - '0+) 

(Mev) R('4+) R('6+) R('8+) 
Po?!0 1.183 1.21 +0.01 1.25 +0.01 spate 
zr*® 2.197 1.41 +0.03 1.57 +0.03 1.64+0.04 
Cr? 1.433 1.65 +0.01 2.17 +0.02 
Sn!” 1.170 1.89 +0.02 2.06 +0.02 _ 
Na‘“4 0.696 1.89 +0.01 2.56 +0.02 or 
Mo* 0.778 2.08 +0.03 3.10 +0.04 — 
Xe!3? 0.534 2.24 +0.02 3.64 +0.03 ica 
Ti*® 0.986 2.33 +0.02 3.39 +0.02 sks 
ca‘? 0.656 2.35 +0.02 3.78 +0.03 
Hg'® 0.41176 2.43 +0.02 3.97 ee 

4.09 +0.03 

Gd'*4 0.1229 3.020 + 0.005 5.83 +0.01 pou 
Ga'5é 0.08897 3.240 + 0.005 6.57 +0.01 re 
Ga'58 0.07956 3. 292 + 0.003 6.774 + 0.006 11.30 + 0.02 
Dy'® 0.087 3.27 +0.01 6.69 +0.01 
Dy'8? 0.0808 3.290 + 0.005 6.80 +0.02 
Eri6é 0.0805 3.290 + 0.005 6.80 +0.08 ae 
Yb'"4 0.0766 3.32 +0.02 6.87 +0.03 eas 
Hf'”6 0.08835 3.290 + 0.005 6.77 +0.13 rom 
Hf!” 0.09319 3.29 +0.01 6.94 +0.03 11.50 + 0.09 
Hf!80 0.09328 3.31 +0.01 6.88 +0.02 11.64 + 0.06 
Os'88 0.13721 3.18 +0.06 6.39 +0.12 me 
Os'% 0.187 2.92 +0.02 5.59 +0.06 8.88 + 0.08 
Ra?* 0.06776 3.10 +0.01 6.14 +0.02 
U2 0.0472 3.31 +0.02 6.80 +0.04 a 
34 0.04350 3.29 +0.02 6.81 +0.03 11.47 + 0.09 
us 0.04528 3.26 +0.03 6.92 +0.12 hat 
Pu238 0.04407 3.310 + 0.005 6.92 +0.02 11.68 + 0.23 
Pu240 0.04287 3.310 + 0.005 6.810 + 0.005 
Cm2#2 0.04212 3.30 +0.01 6.86 +0.12 ‘ 





of Davydov and Filippov as a function of R('4+) 


values seem to be in agreement with experimental 


results. The agreement of the theoretical pre- 
dictions for transition probability ratios and E2, 
M1 admixtures in the 72+ —'2+ transition is quite 
good.*»” It is interesting to note that the smallest 
R('4+) predicted by this theory is 8/3; clearly 
many nuclei have R('4+) smaller than 8/3 (see 
Fig. 1) so that they are not in agreement with 
this theory. 

The “vibrational” nuclei® have average R(*4+) 
value of 2.3° with values ranging from 2.07 to 
2.63 and the closed shell nuclei have R('4+) 
values smaller than 1.9.° The nuclei above the 
“vibrational” nuclei have been called nonaxially- 
symmetric rotational nuclei in Fig. 1. Of all the 
nuclei which are not nonaxially-symmetric rota- 
tional nuclei only Ti** and Sn’ do not fall within 
their experimental error on the empirical curve. 

The classification of nuclei into different types 
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seems not to be essential, at least with respect 
to this system of levels, because all nuclei fol- 
low the same curves which relate the energies 
of these excited states. 

The author wishes to thank F. T. Porter and 
M. S. Freedman for a critical reading of the 
manuscript. 
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ANTIPARTICLES IN SPACE 


Edward W. Hones, Jr. 
Physics Section, Convair, * San Diego, California 
(Received May 11, 1959) 


About a decade ago, Hoyle’ obtained a steady- 
state model for an expanding universe by intro- 
ducing continuous creation of matter into the field 
equations of general relativity. More recently, 
Morrison’ has speculated that the creation of 
matter throughout the universe may take place 
by formation of particle, antiparticle pairs which 
subsequently are segregated by a repulsive force 
between matter and antimatter. Burbidge and 
Hoyle® have estimated that several observed astro- 
physical phenomena (both within our galaxy and 
outside it) might be explained as due to the pres- 
ence of antimatter. But they conclude that the 
ratio of antimatter to ordinary matter in our 
galaxy cannot exceed ~10~’. The purposes of 
this Letter are (a) to point out that a slightly 
more stringent limit for this ratio is suggested 
by cosmic-ray measurements reported by Perlow 
and Kissinger*»® and (b) to observe that relatively 
simple equipment in earth satellites or in space 
probes should serve either to detect antiparticles 
in our galaxy or to define a much more stringent 
limit on their number density than can presently 
be set. 

If one assumes that antiprotons exist through- 
out the galaxy in the limiting number density, 
| n-, of 10~’/cm* estimated by Burbidge and Hoyle 
and if one takes the p, p annihilation cross sec- 
tion used by them and adopts, for protons in the 
galaxy, the number density, n,, of 1 per cm’, 
one obtains a p, p annihilation rate of 3 x10~”” 
em~* sec™. From an average annihilation, there 
ultimately appear 3.6 photons of ~170 Mev and 
3.6 photons of 0.51 Mev.® Therefore, the source 
of gamma rays of each type has an intensity, 
throughout the galaxy, of 10.8x10~*? cm™ sec™. 
The flux, Q (cm~ sec) of these photons coming 
from a spherical region of radius R (cm) sur- 
rounding the earth is 


Q=3PR cm™ sec™, 


. (1) 


where P is the photon generation rate in the 
source (cm™~* sec), We can obtain, simply, a 
very conservative value for the flux, by calcula- 
ting those coming from a sphere of 10000 light 
years (107 cm) radius. (The flux from the whole 
galaxy would certainly be several times larger 
than this.) We find, for each type of photon 

(170 Mev and 0.51 Mev), 


Q = 4(10.8 x 10~*? x 107”) = 5.4 photons/cm? sec. 


In terms of energy flux, this is ~900 Mev/cm? sec 
for the 170-Mev photons and ~2.5 Mev/cm? sec 
for the 0.51-Mev photons. 

Perlow and Kissinger measured the primary 
gamma radiation at 80 km altitude and found the 
flux of photons in the energy range 0.1 to 15 Mev 
to be =2.7/cm? sec from the upper hemisphere.‘ 
This figure seems to be in fair agreement with 
the calculated flux of 0.51-Mev photons until one 
recalls that (a) the calculated flux is not that 
from the whole galaxy but only from a 1078-cm 
sphere and (b) the measurements included not 
only 0.51-Mev photons but all photons between 
0.1 and 15 Mev. From these considerations one 
is probably justified in saying that n-/n,<10-* 
in the galaxy. Perlow and Kissinger measured, 
also, the flux of gamma rays of energy between 
3.4 and 90 Mev and found an energy flux of 1.4 
Mev/cm? sec from the upper hemisphere.® It is 
not possible to make a very meaningful compar - 
ison of these data with the calculated flux of 
170-Mev gamma rays because the latter (if they 
occur) would not have been detected very effi- 
ciently by the apparatus used. However, the 
fact that evidence of this large energy flux (~ 900 
Mev/cm? sec) has not been found in other cos- 
mic-ray experiments indicates that the flux of 
~170-Mev gammas is at least one or two orders 
of magnitude smaller than that calculated above 
(5.4/cm? sec). 
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Considerably more precise information con- 
cerning the density of antiparticles in the galaxy 
can be obtained by carrying a simple scintilla- 
tion gamma-ray spectrometer some distance 
from the earth in a satellite or space probe to 
look specifically for 0.51-Mev photons.’ A search 
could be made, also, for the high-energy (~ 170 
Mev) gamma rays, although the spectrometer to 
be used for these would be larger and more com- 
plicated. 

An explicit search for positrons in space also 
can be made quite simply, using an experimental 
arrangement such as that shown schematically in 
Fig. 1. With this equipment one would detect 
positrons arriving at the annihilation block by 
the resulting coincident 0.51-Mev light pulses in 


e* 


#-— Outer Shell of Satellite 





FIG. 1. 


Schematic representation of satellite instru- 
mentation to detect positrons in space. A-—annihilation 
block of low-Z material; B, C—Nal crystals shielded 
from sun by ~ 1 mg/cm? foil over surface; D, E—lead 
shields to reduce individual counting rates of the two 
Nal crystals; F,G—light pipes; H, J—photomultipliers. 
The system would include electronic circuits which 
would register only coincident 0.51-Mev pulses in 
crystals B and C. 


the two Nal crystals. Using the interplanetary 
cosmic-ray flux found by Van Allen* and assum- 
ing a resolving time of 10~" sec for Nal crystals, 
I estimate that at distances from the earth greater 
than ~ 20000 miles this type of instrument could 
detect positrons if their number density were 
10-**/cm$ or greater. 

Detection of positrons in space in the neighbor- 
hood of the earth’s orbit would not necessarily 
be firm proof of the existence of antimatter on 
a cosmic (or even galactic) scale, because posi- 
trons from cosmic-ray showers or from other 
processes in the solar system might find their 
way to this region. Thus, the search for the 
characteristic annihilation gamma rays probably 
represents the better test for the widespread 
existence of antiparticles. 





"A division of General Dynamics Corporation. 
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K-MESON NUCLEON INTERACTION* 


Robert Karplus, Leroy Kerth, and Thaddeus Kycia 
Department of Physics and Lawrence Radiation Laboratory, | 


University of California, Berkeley, California 
(Received May 15, 1959) 


It is to be expected that the experimental study 
of K-meson production, scattering, and absorp- 
tion processes will shed light on the fundamental 
interactions of the strange particles. The most 
promising method of analysis for determining 
the coupling constants and relative parities of 
the particles exploits the analyticity properties 
of various reaction amplitudes.'»? The recently 
measured angular distribution of the K*-proton 
scattering cross section at various intermediate 
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energies® and the K -proton scattering parame- 
ters*,® enable us to make a more accurate study 
of the K-nucleon-hyperon coupling terms than 
has previously been possible. Unfortunately, the 
accuracy is not sufficient to lead to a definite 
conclusion about the K-meson parity, but we do 
conclude that the K-proton force is attractive. 
We hope that the present discussion also has 
value because it will make clear some of the 
difficulties that must be faced. 








VoL 


Th 
if on 
tion 
= my, 


fle 


wher 
prot 
is th 


Sin 
it is 
does 
used 
the f 


wel 


"> 


in ter 
hyper 

The 
perim 
value 
energ 
nonph 
which 

The 
from | 
the ph 
sectio 




















VOLUME 2, NUMBER 12 


PHYSICAL REVIEW LETTERS 


JuNE 15, 1959 





The numerical work is carried out most easily 
if one uses the observations to calculate the func- 
tion f(w) of K-meson laboratory energy w, (f =c 
=my =1) 


$ 4 
a1 [ A+’) ,,1/ A)... 
ptw)=1[ am *_ e+e dw’, (1) 
1 Wan 
where A,(w) are the imaginary parts of the K*- 


proton forward scattering amplitudes and w,, 
is the (nonphysical) threshold for the reaction 


K°+p-A+n. (2) 


Since the integrals in Eq. (1) have been cut off, 
it is necessary to use a dispersion relation that 
does not weight the high-energy region. We have 
used this expression for the real parts D, (w) of 
the forward scattering amplitudes: 


WyD ,(w) - 3(w, + w)D4(ay) - 3(w - w)D (wy) 
= wpof (w) a (a, + w)f (a) = 3(w, a w)f (-w,) 
PAX 


(w ,” - wy )(w, +w) 


eX > 


2 2 - (3) 
5 7% Mw, + w) 





= Wy (Wy? 9} 





- 


Here p,q (ps) is the parity of K~A (KZ) relative 
to the nucleon and X , (X;), proportional to the 
K~ AN (K~=N) coupling constant g,’ (g5”), is the 
magnitude of the residue at the pole w=w, (ws); 
these quantities are given by the expressions 
(A=A, z) 


” ‘i. «. 2 
w, =(m, *, my )/2m, » 


2 = 
8, (m, +b,m,) my 4) 


40 





4m Pa » 
in terms of the masses M)» Mp; and mr of the 
hyperons, proton, and K meson, respectively. 

The energy w, in Eq. (3) is one at which ex- 
perimental information is available. By using a 
value substantially different from the meson rest 
energy, one can suppress the importance of the 
nonphysical region and of the amplitude D_(w,), 
which decreases with energy. 

The optical theorem is used to obtain A, (w) 
from the measurements of total cross section in 
the physical region. The K*-proton total cross 
sections were taken from the data obtained by 


the use of nuclear emulsions® for energies below 
200 Mev, the results of our experiment,* and the 
data of Burrowes et al.’ for the higher energies. 
A curve was fitted to the data within an estimated 
statistical error of 10%. The low-energy K™- 
proton total cross sections have been measured . 
in a hydrogen bubble chamber.® At higher ener- 
gies the values are o_ =52+ 9 mb at w =2.08,° 
and og. =60+20 mb at w=2.52." We assumed the 
value o. =45 at w =4, and passed a smooth curve 
through all the points. The statistical error of 
the fit was estimated at 15%. 

The values of A_(w) in the nonphysical region 
were obtained by extrapolating with the parame- 
ters of Dalitz and Tuan.**® The integral is very 
sensitive to the sign of D_(1), the real part of 
the scattering amplitude at zero energy.* The 
combination of the integral and the scattering 
amplitude that occurs in Eq. (3), however, is 
not sensitive to this sign as long as the zero- 
range approximation gives the correct sign and 
a reasonably accurate value for D_(w,). This 
is the case because the zero-range formula gives 
an analytic expression for the scattering ampli- 
tude which then essentially satisfies a dispersion 
relation by itself. 

The results for the integration are given in 
Table I. They were obtained with a positive choice 
of D_(1). In the fourth column are the measured 
values® of D,(w); their magnitudes are simply re- 
lated to the total cross sections, because it was 
found that the scattering is isotropic at w =1.46 
(225 Mev) and it was assumed that it is isotropic 
at the lower energies. In the last column is the 
value of D_ obtained from the zero-range param- 
eters, which represent the data adequately.* An 
examination of the bubble chamber results® indi- 
cates that most of the elastic scattering is indeed 
a diffraction effect (i.e., the real part of the 


Table I. Integrals and scattering amplitudes used in 
the dispersion relation Eq. (3). The K-meson energy 
is w measured in units of the rest energy. The unit 
for the other quantities is the K-meson Compton wave- 


length. 











w F() F(-w) D,(w) D_(w) 
1.00 3.0+0.3 -1.25+0.14 
1.17 3.140.3 5.241.0 -1.2420.14 +0.56 
1.285 3.140.3 4.940.8 -1.23+0.14 +0.40 
1.46 3.0+0.3 -1.20+ 0.08 
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scattering amplitude is small) and that there is 
accordingly some tendency toward forward scat- 
tering. Considering the experimental accuracy, 
one may estimate |D_(1.17)|<2 and |D_(1.285)!<2, 
consistent with Table I. 

Since it is clearly not possible to solve for X 
and X5. separately, we have calculated the “aver- 
age” coupling (pX) a, at an “average” pole w 
= 3(w At wy), 


Pax a 


A - w, (w, +w) 





(PX) ay = 3(@ - w,”)(@+ wi 


(w 


PXs 


_ we (w,, =a) (5) 





+ 

(w,, 
for various choices of w and w, in Eq. (3). The 
results are listed in Table II. We conclude that 
this determination of (pX) as leads to the value 


(PX), , =0.0+0.5, (6) 


where we have included some provision for the 
uncertainty in D_(a,). 

Our conclusion, then, is that the K*-proton in- 
teraction, at the present accuracy, is not suffi- 
ciently sensitive to the coupling parameters X A 
and X>5 to permit their evaluation from the pre- 
sent data. This is not too surprising, because 
Eq. (3) depends only on the difference of ampli- 
tudes at closely spaced energies that are quite 
distant from the poles w A and Ws. We must also 
point out an important weakness of this type of 
analysis: only measurements on the K* -neutron 
system can make possible a separation of the 
KAP and K=P coupling terms by their different 
isotopic properties. The determination of the 
K-meson parity from the angular dependence of 
associated KA and KZ photoproduction’ does not 


Table Il. “Average” couplings determined from Eq. 
(3) by the use of data at the energies w and w). The 
results are independent, to the accuracy of this work, 
of the choice of sign for D_(w). 











w Wp PX) ay 
1. 46 1.285 -0.2+0.6 
1.46 1.17 -0.1+0.4 
1.00 1.285 +0.140.4 
1.00 1.17 +0.3+0.6 








512 


suffer from this defect and may, therefore, lead 
to a definitive result sooner.” 

A few words are in order about the amplitude 
D_(w).” If it deviates substantially from its zero- 
range value at the energies of interest, then Eq. 
(6) will be modified somewhat [the right-hand 
side contains D_(w,) with a coefficient -0.2]. The 
characteristic feature of the zero-range approx- 
imation, that the scattering is very strong in 
one of the two isotopic spin states, * depends es- 
sentially only on the fact that the large K~-p 
elastic cross section is accompanied by a sub- 
stantial charge-exchange cross section at low 
energy. Each satisfactory set of parameters in- 
cludes one scattering length whose real part ex- 
ceeds 1.5 fermis in magnitude. We should like 
to add that only an extremely unlikely change in 
the experimental situation could reduce this 
value to 1 f. Since the elastic K~-proton inter - 
action cannot have a range exceeding #/2m_c 
~0.7 f, it follows that the interaction is attrac- 
tive in the isotopic state with the large scattering 
length. If D_(1) is negative, then the attraction 
is so strong as to give rise to a “bound” state 
which will be unstable to decay into a pion-hy- 
peron system. D_(w) can be expected to change 
sign near w~1.2, a meson laboratory momentum 
of 300 Mev/c. If D_(1) is positive, then the 
attraction is not strong enough to give rise toa 
“bound” state. In this case D_(w) will decrease 
monotonically from its low-energy value; near 
w~1.2, it will already be quite small. Since the 
forces are attractive and strong in both of these 
cases, one cannot expect to distinguish between 
them by a study of the nuclear “optical model” 
potential for K~ mesons, which is then not simply 
related to the low-energy scattering amplitude. 
A measurement of the Coulomb interference 
effects in the angular distribution of K -proton 
scattering will, of course, decide the matter. 
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MAGNITUDE OF RENORMALIZATION CONSTANTS* 


S. G. Gasiorowicz 
Lawrence Radiation Laboratory, University of California, Berkeley, California 


D. R. Yennie and H. Suura 
University of Minnesota, Minneapolis, Minnesota 
(Received May 11, 1959) 


Five years ago Kallén' published a proof that 
at least one of the renormalization constants in 
quantum electrodynamics is infinite. A subse- 
quent more refined analysis’ led to the conclusion 
that it is the electron wave function renormaliza- 
tion which diverges. Recently Johnson® raised 
the question whether a formulation of the theory 
in another gauge might not alter this result. We 
have re-examined this problem and arrive at the 
conclusion that (a) Kallén’s proof is not conclu- 
sive, (b) the renormalization constants could be 
finite under rather special circumstances, and 
(c) the question of gauge invariance is quite ir- 
relevant to this problem. We begin with a brief 
review of Kallén’s work on this problem.‘ It can 
be shown that the charge renormalization (1-L)-” 

can be expressed in the form 


(1-L)"*=1+11(0), (1) 





where by definition 


men =P aa Oe, (2) 
and 
V 
11(Q?) = - 30F 2. (O1j,, Im) mij 10); (3) 


P denotes the principal value. Kallén shows that 
in spite of the indefinite metric associated with 
quantum electrodynamics, [I(-a) is positive for 
positive a, and that therefore a lower limit to 
the integral in Eq. (1) can be obtained by con- 
sidering any subset of eigenstates of the Hamil- 
tonian. Thé'simplest of these are the states con- 
sisting of an incoming electron-positron pair. 
Our task therefore is to study the high-energy 
limit of the matrix element (01j,,|p,)’). 

By the use of reduction formulas'’*® one may 
obtain the following compact expression for the 
matrix element in question: 


(01), 1b,b") = (01f, 1p, b’)[1 - (Q*) + (0) - inT1(Q*) + F, (Q*) - F,(0)]+4Q (Olm, 1p,0F,(Q), (4) 


Q =p +p’, 
where 


(4a) 


Mul’ Ifdtxdty ef x4PI (2) 6(-y)X 0 HD), Fe) 1,40} 10) - 6(-»)@(y-x)< 0 {Li , ©), £0), Fle)}10 )}}v@) 


= F,(Q*)(0 i, 10, b") + IF(@E < Olm > P,b’)- (5) 
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In terms of Kallen’s notation, F,(Q*) and F,(Q*) 
are given by® 


F,(Q?) + F,(@*) =F’ 8(Q") + ink” F(@), 
F,(Q@) =5" 8 (Q*) + ins” FQ). (6) 


Here N stands for the electron wave function re- 
normalization constant and f(x) is the fermion 
source density. The matrix elements with the 
superscript (0) are the lowest order perturba- 
tion expressions for the current and the spin. 

In this notation Ward’s identity takes the form 


F, (0) =(1 - M)/(1 - L). (7) 


To estimate the behavior of this matrix ele- 
ment as -Q?~«, Kallén argues as follows: let 
us assume that all the renormalization constants 
are finite. In particular this means that [1(0) is 
finite, so that II(-a) goes to zero fast enough as 
a~«, For any reasonable” function [I(-a), this 
implies that the same must hold for [1(Q*). The 
problem therefore reduces to a study of the 
functions F,(Q*) and F,(Q”). 

Kallén originally examined the left-hand side 
of Eq. (5), and reduced it to an integral of the 





or, more relevantly, 


R™°8(Q2) -R™°F(0) 4 (4Q2)R TE (0)4...4 


da R™ 8-9) 


ere | = 
n 
o +a 


and 5*°8 (Q2) is given by a similar expression. 
The polynomials are seen to be real. It is im- 
portant to note that the degree of the polynomial 
cannot be smaller but may be larger than the 
power necessary to make the integral over the 
imaginary part convergent.® Let us now discuss 
a few possibilities. We shall talk about R but the 
analysis applies equally well to S. 

Case I. n=0.—In this case we have 





_reg eB (_a) 
0 
This implies that 
eo reg 
RTC (gr) RTC Q)--| 222 Ca Gp) 


“a 1+a/@ ? 
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form 


[fo 


As ($,+P,')- this integral need not vanish, be- 
cause p and p’ also vary when p, and p,’ vary 
(since * +m? =p,?; f’?+m?=p,"). The integral 
would vanish if the function F,,@, x; p’, y) ap- 
proached zero as p, p’~« but in general there is 
no way of establishing this. As the vanishing of 
this integral was an important part of Kallén’s 
original proof, we conclude that the proof is in- 
conclusive. 

A more constructive discussion can be given if 
we abandon a certain degree of generality, and 
assume that the functions F,(Q*) and F,(Q”) are 
analytic in the -Q? plane except for a cut along 
the positive real axis.* The behavior of the 
functions F,(@?), F,(Q?) for large values of -Q? 
is unknown. The locality of the theory demands 
that they do not increase faster than a polynomial. 
In general we can therefore write 


1” p,-a) 
F.(Q?) = | 
1 vs 
0 


F > x; p’,y) 
® Teay- Barbs! ~i€) ly -Po’+ie) 





(8) 


@ 


da as" te on (Q?), (9) 


(4?) 


pree- 1) 
LON (0) 


(10) 
a+ 





so that 


lim R™°%(Q*) - RTS (0) - 2 RPE (-g) =0. (13) 


-Q? a 


it 5*°8(Q?) has the same behavior, we obtain 
[using Eq. (7)] the result 


. ’ N ; to) ’ 
(Ol, Ibs") => (014, 19,0"), (14) 


e., the matrix element approaches the unre- 
normalized Born approximation, and a contra- 
diction with the initial assertion that all renor- 
malization constants are finite is obtained be- 
cause [I(-a)—const N*/(1-L)* as a~«. This cor- 
responds to the case discussed by Kallén. 
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Case II. n=1.-—In this case we have 


ig reg 
RPE) «RF 0) -@/ ae ko, (15) 


If the integral 


i <) 


/ ao RTE (a) (16) 
0 


converges, we have 


eo reg 
lim RE) «RE 0) -[ da x (17) 
0 


The difference between this case and the one 
discussed before is that the right side of Eq. (17) 
need not vanish. In particular, if this constant 
is such as to make ( 017, |p, p’) -0 as -Q?—x, 
(this necessitates that there be no polynomial in 
the ST©8 equation), then no conclusion as to the 
magnitude of the renormalization constants can 
be drawn from the electron-positron pair states 
alone. This is the only possible exception to 
Kallén’s result. If the integral in Eq. (16) does 
not converge, R*™°8(Q?) will diverge in the limit 
-Q?~o for any reasonable function R™°8(-a). 
Thus except for the special case noted above, 
the integral in Eq. (2) will always diverge (this 
is even more easily seen if n=2,3,...).'° The 
special case corresponds to the “no-subtrac- 
tion” spectral representations for the matrix 
elements of the current operator discussed 
recently." 

It is clear from this discussion that developing 
the theory in another gauge cannot make any 
difference. We have used gauge invariance only 
in the particular form in Eq. (5), but the scalar 
functions F,(Q?) and F,(Q?) can be discussed in- 
dependently of whether the photon has a mass or 
not. This fact will affect the low-frequency be- 
havior of the functions [e.g., Eq. (7)] but cannot 
make any difference in the high-energy region. 
Our negative conclusion applies by the same 
token to all renormalizable theories: one cannot 
prove that the renormalization constants are in- 
finite, because there is always one possibility 
that a special cancellation occurs. 

We may ask whether it is possible that the re- 
normalization constants are indeed finite. For 
this to be true all matrix elements in Eq. (3) 
must have the “no-subtraction” form. Actually 
this is not enough, because even though every 
contribution to the integral [1(0) will be finite, it 





is not necessarily true that the sum of these con- 
tributions converges. In fact, since (0 ly |n) could 
not be expected to approach zero until a>>(E,N)? 
(where E, is some fixed energy and N is the 
number of particles in the state), a very nonuni- 
form series is involved, and counterexamples to’ 
convergence can readily be constructed.** Thus 
unless field theory has some regularity proper- 
ties quite unsuspected from perturbation theory, 
the renormalization constants will be infinite. 

In our discussion we have given up a certain 
amount of generality by working with the un- 
proved spectral representations. Should these 
turn out to be wrong, there may be two possibil- 
ities: If the violation of analyticity is limited to 
a finite part of the plane (e.g., longer cuts), then 
the domination of the high-energy behavior by a 
polynomial is unaltered and our discussion in 
terms of spectral representations needs only 
minor modifications. This would not be true, 
however, if regions of nonanalyticity extend to 
infinity (e.g., cuts not indicated by perturbation 
theory such as those which might arise from 
complex “ghosts”). 

This work was started while one of us (S.G.) 
was visiting the Department of Physics at the 
University of Minnesota. He is grateful to Pro- 
fessor A. O. Nier and the Department of Physics 
for their kind hospitality. 
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PRECISE DETERMINATION OF THE MUON 
MAGNETIC MOMENT. R. L. Garwin, D. P. 
Hutchinson, S. Penman, and G. Shapiro [Phys. 
Rev. Letters 2, 213 (1959) ]. 


It has recently been’ established that the theo- 
retical uncertainties in the lower limit of the 
muon mass from mesonic a-rays? are much 
smaller than previously supposed.* We there- 
fore undertook a more careful search into pos- 
sible systematic errors in our moment measure- 
ment. We have found and corrected one serious 
source of error. A method of testing was then 
designed which simulated the experimental situa- 
tion and permitted the establishment of an upper 
limit on the remaining systematic errors in the 
electronic equipment of 1 x10~* or almost an 
order of magnitude less than our stated accuracy. 

The circuit design used in this experiment has 
been intended to be completely aperiodic, i.e., 
should introduce no systematic errors in a fre- 
quency measurement (as opposed to a random 
error in measuring each particle due to the fi- 
nite time resolution). Preceding the analysis 
equipment, however, is the “zero-crossing de- 
tector” which does the fast timing of both the 
muon and decay electron pulses. To circumvent 
the difficulties of long-term transit-time varia- 
tions in phototubes, a single counter and timing 
circuit was used for both muon and electron 
pulses. This raises the possibility that the re- 
covery of the phototube or fast timing circuit 
after the passage of a muon pulse will alter the 
apparent time of the electron. This effect is im- 
portant only if it varies during the measurement 
interval. In particular, if electrons immediately 
following the muon pulse are delayed more than 
those coming later in the gate interval, the net 
effect is to make the apparent muon precession 
frequency appear larger than it actually is. Un- 
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fortunately, we observed on our circuit diagram 
a one-microsecond time-constant at the grid 
following the distributed zero-crossing ampli- 
fier. Had this time-constant been 50 millimicro- 
seconds or >20 microseconds, no error would 
have been incurred, but we estimated the maxi- 
mum error caused by it as ~4x10"*. 

Since there is available to us no measuring 
equipment that approaches the experimental ap- 
paratus in resolving ability, the existence of the 
postulated effect was established by inserting 
two pulses with fixed cable delay of 2 micro- 
seconds into the equipment. The first was ana- 
lyzed by the circuitry as a muon pulse and the 
other as an electron pulse. The equipment out- 
put was displayed on a pulse-height analyzer as 
previously described. A third pulse which had 
been traveling in a long cable was then inserted 
between the “muon” and “electron” pulses. Its 
retarding effect on the electron timing could then 
be seen clearly. The fast timing circuit was then 
altered (by changing the offending time-constant) 
so that the effect was no longer apparent. 

After modifying the circuitry, the following 
rigorous test was made to establish that there 
remained no other such sources of error. We 
built a source of random pulses whose probabil- 
ity of occurrence above a fixed threshold oscil- 
lated in time at a known rate. This consisted of 
a 6810A photomultiplier viewing a plastic scintil- 
lator which was irradiated with 8 rays. The 
focusing voltage of the tube was modulated at or 
near the frequency of the reference oscillator of 
the measuring apparatus, (86.2 Mc/sec). A 30- 
volt peak-to-peak signal was sufficient to a- 
chieve nearly 100% modulation of the output pul- 
ses. The slow coincidence circuits select the 
first pulse above a threshold after a dead time 
as a yu pulse and the next puise as an electron. 
Since the probability of occurrence of these two 
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pulses is strongly correlated at the reference 
oscillator frequency, the pulse-height analyzer 
display of the electron phase modulo the refer- 
ence oscillator period will show a structure very 
similar to that for a precessing muon [although 
not restricted to the form (1+a cos@)]. The ac- 
tual experiment consists of comparing the phase 
distribution of electrons that come soon after a 
muon stop with those that come later. In the 
test, two groups of modulated random pulses 
were displayed, those that corresponded to “early” 
electrons and those that corresponded to “late” 
ones. 

Figure 1 shows the actual pulse-height analyzer 
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FIG. 1. Pulse-height analyzer display of the mod- 


ulated counter test. 


display. On the left is the distribution in phase 
of the early “electrons” and on the right that of 
the late “electrons.” Ideally, the distributions 
will be identical in shape. A treatment of the 
test data by the same Fourier analysis used in 
the experiment shows no phase shift correspond- 
ing to a frequency error larger than 10°°. 

The measurement of the muon moment was then 
repeated using a smaller bromoform target, with 
the equipment modified and tested as above. 

Expressing the final result in terms of the 
ratio of the resonance frequency of the muon to 
that of a proton in the same magnetic field, we 
obtained f,, /fp = 3.1834 +0.0002. Combining this 
value for the moment with the value for the 
lower limit on the mass, as discussed in refer- 
ence 1, we obtain for the g factor g,, >2(1.00122 
+0.00008). This limit is seen to be in agreement 
with the theoretical prediction g=2(1.00116). The 
experiment to determine the absorption coeffi- 
cient of lead for the 3D - 2P transition in u-mes- 
onic phosphorus is being repeated. Parratt* has 
pointed out that the width of the lead K edge 
should enable one to obtain either a more accu- 
rate mass value or a considerably different 
lower limit. 





1A, Petermann and Y. Yamaguchi, Phys. Rev. Let- 
ters 2, 359 (1959). 

*Koslov, Fitch, and Rainwater, Phys. Rev. 95, 291 
(1954). 

8Cohen, Crowe, and DuMond, Fundamental Constants 
of Physics (Interscience Publishers, New York, 1957). 

‘L. G. Parratt (private communication). 
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ABSTRACTS meaningful even for potentials with hard cores, a low 
and a nonzero solution implies the existence of pecte 
In this section are printed the abstracts of Articles that a superfluid state. A variational method is de- of the 
have been forwarded to The American Institute of Physics vised which provides a criterion for superfluidity By 
for publication in THE PHYSICAL REVIEW’. In quoting and a lower bound for the transition temperature | face- 
information obtained from this section before the appear- into the superfluid state. We find that a repul- beam 
ance of the corresponding Article, reference should he sive hard core does not in principle forbid the be ob 
made to “Physical Review (to be published)” rather than existence of a superfluid state, but whereas in of the 
to this Journal. the absence of a hard core an attractive two- thern 
body potential always leads to a superfluid state inter« 
MEASUREMENT OF DISTANCE IN GENERAL at sufficiently low temperatures, in the presence 
RELATIVITY. E. Newman, University of Pitts- of a repulsive core there appears to be a critical 
burgh, Pittsburgh, Pennsylvania, and J. N. Gold- strength of attraction needed to form a super- OPTI 
berg, Aeronautical Research Laboratory, Wright- fluid state. When the variational principle is ap- Newn 
Patterson Air Force Base, Ohio (Received plied to liquid He® or to nuclear matter, it is Rese: 
January 28, 1959). found for a wide class of trial functions that the (Rece 
An intrinsic geometrical definition of distance SPEER SONS ES HESSEN & GNSS. The 
is presented in terms of the geodesic deviation “Present address: Brown University, Providence, has b 
of null rays. This definition is applied to the Rhode Island. its re 
Schwarzschild solution and the homogeneous and The z 
isotropic cosmological solutions of the Einstein STATIONARY NONEQUILIBRIUM GIBBSIAN from 
field equations. In the former case the distance ENSEMBLES. J. L. Lebowitz, Department of of the 
is proportional to the radial coordinate of the Physics, Stevens Institute of Technology, Hobo- tion ¢ 
standard metric, and in the latter case it yields ken, New Jersey (Received January 6, 1959). magn 
the usual cosmological distance. However, since The : 
the geodesic deviation is a scalar, the results are The general theory of a Gibbs ensemble repre- ]| g « 
independent of the particular metric used. The senting a system in contact with its surroundings | .1.., 
relationship of this definition to observation is is applied to several concrete situations of in- Sound 
discussed and it is concluded that it agrees with terest. By an appropriate choice of heat reser - sity « 
the astronomical definition of luminosity distance. VOirs a simply modified Liouville equation is behat 
found to describe a heat conducting system. The antife 
stationary nonequilibrium I'-space ensembles spect 
which describe such a system are found explicitly eters 
POSSIBLE SUPERFLUIDITY OF A SYSTEM OF for some cases. In the simplest cases these iain 
STRONGLY INTERACTING FERMIONS. L. N. ensembles turn out to be canonical with a temper- 7 ve 
Cooper,” R. L. Mills, and A. M. Sessler, Ohio ature that is a weighted average of the reservoir a 
State University, Columbus, Ohio (Received temperatures. For other systems, such as ome 
January 30, 1959). Brownian particles inside a fluid whose temper - 
, —_ ature is not uniform, we find the stationary en- STRE 
The possible superfluidity of a system of Bore THE 
strongly interacting fermions is investigated on semble to terms linear in the temperature gra- OF N 
the assumption that an adequate description of dient. Fr _ this ~~ e led to discuss — Seals 
the system in its “normal” state is given by in- bles that will approximately represent an arbi- Reads 
dependent fermions in a momentum -dependent trary heat conducting fluid. A more general 1959) 
potential. On the basis of this assumption we proof than previously given for the asymptotic 
have investigated whether a correlated wave func-  2PProach of the I'-space distribution to its sta- The 
tion of the form used by Bardeen, Cooper, and tionary value is also presented. state 
oe minimizes the ground-state energy. VELOCITIES OF SPUTTERED ATOMS. G. K. ei 
e nonzero terms in the expectation value of , a , Pp 
2 : . ae Fg Wehner, Mechanical Division of General Mills, 
the Hamiltonian contain the modified kinetic en- , : ; to 10! 
é Incorporated, Minneapolis, Minnesota (Received 
ergy and the full two-body potential between the iunuary 6, 1000) The a 
fermion pairs. An integral equation is obtained ’ . pansir 
in configuration space for the correlation func - The average velocity of atoms sputtered from termi 
metal surfaces by normally incident Hg* ions in know! 


tion between pairs. This integral equation is 
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a low-pressure discharge was found to be unex- 
pectedly high, corresponding to kinetic energies 
of the order of 10 to 30 ev. 

By sputtering atoms from a (110) surface of a 
face-centered cubic crystal, rather dense atomic 
beams ejected normal to the target surface can 
be obtained. The fact that the kinetic energies 
of the atoms are more than 100 times higher than 
thermal evaporation energies makes such a beam 
interesting. 


OPTICAL PROPERTIES OF NICKEL OXIDE. R. 
Newman and R. M. Chrenko, General Electric 
Research Laboratory, Schenectady, New York 
(Received January 28, 1959). 


The absorption spectrum of single-crystal NiO 
has been measured in the 0.1 to 6 ev range and 
its reflectivity spectrum from 0.025 to 10 ev. 

The absorption spectrum shows a series of lines 
from 1 to 3.5 ev arising from internal transitions 
of the Ni ion. A continuous background absorp- 
tion occurs in the range from 0.1 to 3.5 ev, whose 
magnitude increases with impurity concentration. 
The absorption coefficient rises steeply above 

3.5 ev and reaches a value of 10° cm~ at and 
above 4 ev. An absorption line at 0.24 ev is 
found to be temperature sensitive in both inten- 
sity and frequency in the range above 300°K. Its 
behavior suggests that it is connected with the 
antiferromagnetic ordering. The reststrahlen 
spectrum was observed with the following param - 
eters: high- and low-frequency dielectric con- 
stants 5.4 and 12, respectively; energies of 
longitudinal and transverse optical mode vibra- 
tions are 0.076 and 0.044 ev, respectively. 


STRESS AND TEMPERATURE DEPENDENCE OF 
THE PARAMAGNETIC RESONANCE SPECTRUM 
OF NICKEL FLUOSILICATE. W. M. Walsh, Jr., 
Gordon McKay Laboratory, Harvard University, 
Cambridge, Massachusetts (Received January 19, 
1959). 


The crystalline field splitting of the ground- 
state spin triplet of divalent nickel in the fluosili- 
cate, NiSiF,-6H,O, has been measured at room 
temperature as a function of hydrostatic pressure 


to 10000 kg/cm? and uniaxial stress to 150 kg/cm’. 


The anisotropic compressibility and thermal ex- 

pansion of this trigonal crystal have also been de- 
termined. When these data are combined with the 
known variation of the splitting with temperature, 


the dependence of the splitting on isothermal unit- 
cell geometry and on temperature at constant unit- 
cell dimensions is calculated. The splitting is 
found to be independent of volume within experi- 
mental error but proves to be quite sensitive to 
unit-cell shape. The deduced explicit tempera- - 
ture dependence is three times larger than that 
measured at atmospheric pressure. The magni- 
tude and geometrical variation of the crystalline 
field splitting may be qualitatively understood 
using a static, ionic model of the (Ni- 6H,O)?* 
octahedral complex. A rather general analysis 

of the explicit temperature dependence indicates, 
however, that low-frequency lattice vibrations 
play a dominant role in determining the observed 
value of the splitting. 

The resonance line widths are observed to in- 
crease monotonically and quite nonlinearly with 
increasing pressure. This broadening is discussed 
in terms of isotropic and anisotropic exchange 
interactions. In agreement with earlier conclu- 
sions of Ollom and Van Vleck, it is inferred that 
the two mechanisms are of comparable impor - 
tance in this paramagnetic salt. 


PRESSURE DEPENDENCE OF THE PARAMAG- 
NETIC RESONANCE SPECTRA OF TWO DILUTE 
CHROMIUM SALTS. Walter M. Walsh, Jr., 
Gordon McKay Laboratory, Harvard University, 
Cambridge, Massachusetts (Received January 19, 
1959). 


The pressure dependences of the spin-Hamil- 
tonian parameters of trivalent chromium in am- 
monium aluminum alum and potassium cobalti- 
cyanide have been measured up to 10000 kg/cm? 
near room temperature. In the case of the alum, 
runs have been made at 0, 24, and 50°C. The 
g-value remains unchanged while the crystalline 
field splitting, 5, increases by ~30%, (85/aP)7 
decreasing with rising temperature and pressure. 
The compressibility of this cubic crystal is known 
but it was necessary to measure the thermal ex- 
pansion coefficient in order to complete the equa- 
tion of state. The crystal contracts when heated: 
(alnV /8T) = -8.1 x 1075(°C)™ in the interval (0, 50°C). 
This anomalous result is consistent with the tem- 
perature dependence of the compressibility found 
by Bridgman and may correlate with the known 
low-temperature phase transition of the alum. 
Using the empirically determined equation of state, 
the crystalline field splitting variations are con- 
verted to isothermal volume and explicit tem - 
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perature dependences. These are discussed in 
terms of static and dynamic crystalline fields, 
but are not satisfactorily explained. This failure 
is attributed to inhomogeneous internal deforma- 
tion of the unit cell as a function of stress and 
temperature. 

A room temperature run on the covalent chromi- 
cyanide shows the g-value as well as the prin- 
cipal splitting parameter, D, to be nonmonotonic 
functions of pressure. The rhombic splitting 
parameter, E, increases quadratically with pres- 
sure. Since no attempt was made to determine 
the crystalline equation of state, the resonance 
data are only qualitatively discussed. It is dif- 
ficult to reconcile the results with the equivalent 
crystalline-field model of the chromicyanide 
complex. 


THERMOLUMINESCENCE AND F CENTERS. I. 
THEORY. G. Bonfiglioli, P. Brovetto, and C. 
Cortese, Istituto Elettrotecnico Nazionale, Turin, 
Italy (Received December 22, 1959). 


After a critical review of the present state of 
knowledge about thermoluminescence, which is 
shown to be rather unsatisfactory, the present 
paper considers in detail thermoluminescence 
arising from thermal bleaching of F centers in 
irradiated alkali halides. A model is proposed, 
which leads to second-order kinetics for the 
radiative recombination of F centers and Lumi- 
nophor centers (L;) which must be present in 
the crystal lattice. 

The formalism concerned is developed, show- 
ing how activation energy and probability factors 
can be determined. 

A quantitative argument is presented, con- 
cerned with the nature of L. centers and suitable 
for the identification of their ground levels with- 
in the forbidden crystal band. 


THERMOLUMINESCENCE AND F CENTERS. IL. 
EXPERIMENTAL. G. Bonfiglioli, P. Brovetto, 
and C. Cortese, Istituto Elettrotecnico Nazionale, 
Turin, Italy (Received December 22, 1959). 


The theoretical treatment developed in the pre- 
vious paper is applied to interpret experiments 
of thermoluminescence associated with the pres- 
ence of F centers in irradiated alkali halides. 
Some important features of the technique em- 
ployed are discussed—in particular, the pro- 
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cedure used to establish the proper wavelength 
belonging to the light emitted in each glow peak. 

The energy of activation and the probability 
factors of the phenomenon are determined, and it 
is shown that the second and third peaks of 
thermoluminescence are probably due to direct 
radiative recombination of F centers with V, and 
V, centers. On the other hand, it is still uncer- 
tain whether the first peak is associated with V, 
centers or with Cu ion impurities in the NaCl 
lattice. 


TEMPERATURE DEPENDENCE OF THE ELEC- 
TRICAL PROPERTIES OF BISMUTH-ANTIMONY 
ALLOYS. A. L. Jain, Department of Physics and 
Institute for the Study of Metals, University of 
Chicago, Chicago, Dlinois (Received January 22, 
1959). 


The electrical resistivity p and Hall effect R 
of zone-levelled single crystals of Bi-Sb alloys 
have been measured in the temperature range 
from 4.2°K to 300°K. Log(p/p,..) vs 1/T curves 
suggest thermal activation of carriers in the 
concentration range from 5 % to 40% Sb in the 
temperature range from 25°K to 100°K; approxi- 
mate activation energies have been inferred 
from their slopes. The activation energy ap- 
pears to have a maximum at a concentration néar 
12%. Some anomalies have been observed in the 
behavior of p and R on both sides of this concen- 
tration at low temperatures. Lattice parameters 
for these alloys have also been measured for the 
entire range of solid solubility. Both a maximum 
and minimum in the c-axis lattice parameter vs 
concentration occur near the concentration at 
which anomalies appear in transport properties. 
These phenomena are discussed in terms of a 


simple band model proposed by Blount and Cohen. 


ELECTRON BOMBARDMENT OF SILICON. D. E. 


Hill,* Purdue University, West Lafayette, India- 
na (Received January 22, 1959). 


The electrical properties of single-crystal sili- 
con have been investigated as a function of bom- 
bardment with high-energy electrons. Measure- 
ments of the initial carrier removal rate were 
carried out at 273°K for samples with widely 
different carrier concentrations. A second type 
of experiment involved the measurement of the 
temperature dependence of resistivity and Hall 
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coefficient at different stages of the bombard- 
ment. Bombardment-produced energy levels 
were found at 0.03 ev, 0.17 ev, and 0.4 ev below 
the edge of the conduction band. With 4.5-Mev 
electrons, the introduction rates for these levels 
were 11 cm™, 0.5 cm™, and 0.05 cm~, res- 
pectively. Bombardment levels were also found 
at 0.05 ev and 0.3 ev above the edge of the va- 
lence band. With 4.5-Mev electrons, the intro- 
duction rates for these levels were 13 cm™ and 
roughly 0.3 cm™, respectively. Irradiations 
with 700-kev electrons indicate the presence of 
the same defect levels with correspondingly 
smaller introduction rates. An investigation of 
the irradiation-induced mobility changes for de- 
generate p-type silicon indicates that, for each 
hole removed from the valence band, there ap- 
pears a singly-charged scattering center. For 
near-degenerate n-type samples, the mobility 
changes are much smaller. Results of bombard- 
ments at 78°K and 10°K show that the net damage 
rate decreases by a factor of 2 or more when 
bombardments are carried out at lower tempera- 
tures. 

*Now at Monsanto Chemical Company, Dayton, Ohio. 


CYCLOTRON ABSORPTION IN METALLIC BIS- 
MUTH AND ITS ALLOYS. J. K. Galt, W. A. Yager, 
F. R. Merritt, and B. B. Cetlin, Bell Telephone 
Laboratories, Murray Hill, New Jersey, and 

A. D. Brailsford, Ford Motor Company, Dearborn, 
Michigan (Received January 21, 1959). 


Extensive studies of cyclotron absorption in 
bismuth and in alloys of bismuth with tin and tel- 
lurium are reported. Experimentally, a micro- 
wave cavity arrangement is used which causes 
circularly polarized radiation to be incident on 
the sample. The samples are disks which form 
part of the cavity wall during the experiment. A 
signal is observed which is under most condi- 
tions proportional to the change in the fractional 
power absorption of the sample as the magnetic 
field applied to the sample is changed. The ex- 
periment is done with this field in the plane of 
the sample as well as normal to it. When the 
magnetic field is normal to the sample surface, 
circularly polarized radiation distinguishes bet - 
ween the effects of holes and of electrons. Making 
this distinction is complicated in bismuth by the 
anisotropy of the band structure, but it is not 
impossible. In pure bismuth the number of holes 
equals the number of electrons; alloying with tin 


increases the relative number of holes and al- 
loying with tellurium increases the relative num- 
ber of electrons. The experiments have been 
done both at 24000 Mc/sec and at 72000 Mc/sec. 
The tilted ellipsoid model of the electron band, ” 
and the ellipsoid of revolution model for the hole 
band®* have been used in a classical magneto -ionic 
theory of the effect which fits the data taken with 
magnetic field normal to the sample surface quite 
satisfactorily. The fit between this theory and the 
experiment is substantial confirmation of the 
validity of this band model, and makes it pos- 
sible with the aid of certain infrared observations 
to determine the parameters of the model numer - 
ically. The mass parameters obtained for an el- 
lipsoid in the electron band are m, = 0.0088, m, 

= 1.80, m, = 0.02323, m,=+0.16, while those for 
the hole band are M, = M, = 0.068 and M, = 0.92. 
Comparison is made with similar results obtained 
from magneto-resistance and de Haas - van Alphen 
experiments. The experiments done with the mag- 
netic field in the plane of the sample have been 
used to improve the accuracy of the determination 
of cyclotron resonance fields and therefore of the 
masses of the charge carriers, although no quan- 
titative theory of the over-all shape of the signal 
under these conditions is available. 


NUCLEAR MAGNETIC RESONANCE ACOUSTIC 
ABSORPTION IN KI AND KBr. D. I. Bolef and 
M. Menes, Physics Department, Westinghouse 
Research Laboratories, Beulah Road, Churchill 
Borough, Pittsburgh, Pennsylvania (Received 
January 28, 1959). 


A technique for measuring the coupling between 
nuclear spins and the crystal lattice is described. 
The technique is one in which the absorption of 
energy by the nuclear spins from an acoustic 
wave is observed as an additional attenuation of 
the wave in the sample. This attenuation results 
in a change in the electric impedance of the 
piezoelectric transducer which drives the sam- 
ple. The transducer impedance is made to con- 
trol a Pound-Watkins type of spectrometer 
through an inverting and matching network. 

Nuclear magnetic resonance acoustic absorp- 
tion measurements were made on single crystal 
samples of KI and KBr using compressional 
waves in the [100] direction. These measure- 
ments yield a value of the electric quadrupole 
coupling e?4,Q, where q, is the first derivative 
with respect to strain of the electric field gradi- 
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ent at the nuclear position. For I'*’ in KI, 
e74¢,Q/h=300 Mc/sec. For Br™ in KBr, e7¢,Q/h 
= 82 Mc/sec. Interpreting these results on the 
basis of a point charge model, we obtain values 
of the amplification factor y. For KI'*’, y=38. 
For KBr®, » =26. 


MILLIMETER-WAVE ROTATIONAL SPECTRUM 
OF NO IN THE *z,, STATE. Paolo G. Favero,* 
Anna M. Mirri, and Walter Gordy, Department 
of Physics, Duke University, Durham, North 
Carolina (Received January 26, 1959). 


The J=3/2~5/2 rotational transition of NO in 
the *n,,. electronic state has been measured in 
the 1.17-mm wave region. Theory applied to 
these measurements combined with previous 
measurements on transitions of the ground 77, 
state lead to the following values for the charac- 
teristic constants of N**O"*: The rotational con- 
stant, B,=50 848.42 Mc/sec, and the spin-orbit 
coupling constant A =122.094 cm™*. The N** nu- 
clear magnetic couplings a= 83.82 Mc/sec, b 
= 68.49 Mc/sec, c = -86.34 Mc/sec, and the N* 
nuclear quadrupole couplings eQq, = -2 Mc/sec, 
e€Qq,=22 Mc/sec. The A-doubling constants p A 
= 175.15 Mc/sec, GQ, =1-15 Mc/sec. 


* 
Permanent address: Institute of Physical Chemistry, 
University of Padua, Padua, Italy. 


SCATTERING OF ELECTRONS WITH AN INTRIN- 
SIC ELECTRIC DIPOLE MOMENT. B. Margolis, 
S. Rosendorff, and A. Sirlin, Department of Phys- 
ics, Columbia University, New York, New York 
(Received January 19, 1959). 


The elastic scattering of electrons with a small 
intrinsic electric dipole moment by spinless tar- 
gets is discussed in the first and second Born 
approximations. For large momentum transfer 
(¢ ~500 Mev/c) the electric dipole moment con- 
tribution to the elastic scattering cross section 
becomes important even for values 4 ~10~° (where 
4 is the electric dipole moment in units of e times 
the Compton wavelength of the electron). Thus, 
it seems likely that experiments involving high- 
energy electrons may give valuable information 
on the hypothetical existence and upper limit of 
the electric dipole moment interaction. It is also 
found, in the first Born approximation, that at 
sufficiently high energies the existence of the 
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electric dipole moment can change the polari- 
zation of the electron beam from longitudinal to 
transverse. Finally, asymmetry effects which 
appear in the second Born approximation are dis- 
cussed in detail. Detection of these asymmetries 
or of the spin correlation effects of the first Born 
approximation mentioned above would imply the 
noninvariance of the electromagnetic interactions 
under parity and time reversal. 


IONIZATION PRODUCED BY ATOMIC COLLI- 
SIONS AT kev ENERGIES. II. A. Russek and 

M. Tom Thomas, t The University of Connecticut, 
Storrs, Connecticut (Received January 26, 1959). 


The electron evaporation model, presented in 
a previous paper, of the collision-ionization 
process that occurs when atoms collide at high 
energies is extended and improved. In addition, 
an alternate model of the process is considered. 
This second model treats the collision-ioniza- 
tion process as being due to direct knock-outs of 
the electrons by violent electron-electron colli- 
sions. It is found that this model cannot be made 
to account for the data, thereby lending additional 
support to the assumptions inherent in the evapor- 
ation model. 


tPresent address: Physics Department, Brown 
University, Providence, Rhode Island. 


DECAY OF Zn™. James B. Cumming, Chemistry 
Department, Brookhaven National Laboratory, 
Upton, New York (Received January 26, 1959). 


The 88.5-sec activity observed among the pro- 
ducts of the bombardment of nickel with 16-Mev 
alpha particles is assigned to Zn™. The decay 
properties of this isotope have been investigated 
by beta and gamma scintillation techniques and by 
coincidence measurements. Approximately 80% 
of the Zn® beta decays are to the ground state of 
Cu”. The positron end point is 4.38 Mev which 
corresponds to a log/ft of 5.2 and to an allowed 
transition. Gamma rays of 0.48, 0.98, and 1.64 
Mev have been observed with intensities of 11.2, 
2.9, and 6.2 percent of the positrons. There is 
some indication of gamma rays at 0.69 Mev and 
possibly other energies in lower intensities. A 
decay scheme is proposed which incorporates the 
above data into an existing level scheme of Cu™ 
derived from (p,)) experiments. 
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DECAY OF Co and Mn®. P. Kienle,* Brook- 
haven National Laboratory, Upton, New York, 
and R. E. Segel, Aeronautical Research Labora- 
tory, Dayton, Ohio (Received January 23, 1959). 


Scintillation spectroscopy and coincidence tech- 
niques have been used to measure the decay 
schemes of Mn*® and Co, both of which decay to 
Fe*®. Gamma-ray energies and intensities were 
measured with the aid of a large Nal(T1) crystal 
and a three-crystal pair spectrometer. Combin- 
ing these data with y-y and B-y coincidence meas- 
urements, decay schemes are constructed for 
both nuclei. The logft values for the various 8 
branches are examined and all of the 8 branches 
from Co™ show higher than usual logft values 
for allowed transitions. The Co decay to the 
second excited state in Fe® is inhibited by a 
factor of ~10° while the Mn*® decay to the same 
state is inhibited even further. The allowed Mn*™® 
decay to the first excited state in Fe® is also 
slowed, while Mn®® decays to higher states pro- 
ceed at more normal speeds. A discussion of 
the 8-transition rates in terms of shell model 
assignments is given, from which some explana- 
tion of the 8 rates results. The high inhibitions 
of the decays to the second excited state in Fe™® 
do not seem susceptible to a simple shell model 
explanation. 


"Permanent address: Technische Hochschule, Mun- 
chen, Germany. 


THEORY OF DIRECT-INTERACTION INELASTIC 
SCATTERING. Norman K. Glendenning, Lawrence 
Radiation Laboratory, University of California, 


Berkeley, California (Received January 12, 1959). 


A model is proposed for the description of di- 
rect-interaction inelastic scattering in which re- 
solved final states of the target nucleus are ob- 
served. This model is expected to be useful in 
deducing the spins and parities of the excited 
states of nuclei. 


INTERNAL CONVERSION IN Li’. R. D. Leamer,* 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received January 26, 1959). 


Upper and lower limits have been placed on the 
internal conversion coefficient of the 479-kev 
excited state of Li’. The limits are 16.6x10~ 
and 5.8x10~’, yielding a multipole mixture of 


(100% - 55%) M1 and (0% - 45%) £2 for the gamma 
ray. Photoelectron emission from the source and 
the background counting rate impose the princi- 
pal limitations on the accuracy of the experiment. 


“Now at Westinghouse Reactor Evaluation Center, 
Ruffsdale, Pennsylvania. 


FINITE-RANGE DELTA-FUNCTION POTENTIAL. 
M. Bolsterli, School of Physics, University of 
Minnesota, Minneapolis, Minnesota (Received 
January 26, 1959). 


The long-range part of the nucleon- uucleon 
potential is approximated by the finite-range 
delta-function potential: 


Vsinglet even =V,5(r - a), 
Vtriplet even =V35(r -a)+V75(r -a)S,,. 


The best fit to the low-energy n-p data is ob- 
tained when V, is repulsive and V7 is large. The 
uncertainties in the V’s are of the order of 10%. 


NEUTRON ACTIVATION AT 195 kev. W. S. 
Lyon and R. L. Macklin, Oak Ridge National 
Laboratory, Oak Ridge, Tennessee (Received 
January 19, 1959). 


Neutron activation cross sections have been 
measured at 195 kev (50 kev energy spread) for 
about 30 isotopes, by absolute gamma counting. 
Absolute calibration was effected by determining 
the RdTh-D,O neutron source strength with a 
flat-response 47 graphite-sphere neutron counter, 
and irradiating indium in the form of a thin, un- 
iform, spherical shell. The results are in gen- 
eral agreement with those of other investigators. 
The 25-kev activation cross sections of Pr, 
Zr™, Ga”, and Pd’ have been revised. 


EXPERIMENTAL DETERMINATION OF THE 
NONMAGNETIC NEUTRON-ELECTRON INTER- 
ACTION. E. Melkonian, B.M. Rustad, and W. W. 
Havens, Jr., Columbia University, New York, 
New York, and Brookhaven National Laboratory, 
Upton, New York (Received January 21, 1959). 


A precise set of measurements of the variation 
of the transmission of liquid bismuth with energy 
in the range 0.1 ev to 10 ev has been made. After 
correction for neutron capture, Doppler effect, 
and liquid diffraction, a value of V, = -4.34+0.14 
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kev is obtained as the depth of the potential well 
having a radius equal to the classical electron 
radius. This gives -0.26+0.14 kev as the strength 
of the intrinsic neutron-electron interaction after 
subtraction of the Foldy term. 


SELF-INDICATION RESONANCE NEUTRON EX- 
PERIMENTS WITH SILVER (2-200 ev). James E. 
Draper, * Brookhaven National Laboratory, Upton, 
New York, and Yale University, New Haven, 
Connecticut (Received June 5, 1958; revised 
manuscript received January 19, 1959). 


A method is given for the area analysis of res- 
onance neutron self-indication experiments inclu- 
ding consideration of the effects of energy resolu- 
tion, wing area, energy dependence of neutron 
flux, v~ term in the capture cross section, poten- 
tial scattering, detailed treatment of Doppler 
broadening, and scattering events preceding cap- 
ture in the detector sample. The essential result 
of this measurement alone is nearly the peak total 
cross section of the Doppler -broadened resonance 
(viz., O,9¥%9/~9"*) which, with approximate know- 
ledge of [', gives 049. Experimental results are 
given for eleven resonances below 200 ev in sil- 
ver, and they are compared with results obtained 
by other methods. 


- 
Present address: Yale University, New Haven, 
Connecticut. 


METHOD FOR THE DETERMINATION OF 
STRANGE PARTICLE PARITIES AND COUPLING 
CONSTANTS. John G. Taylor, Christ’s College, 
Cambridge, England (Received January 7, 1959). 


A general method for determining parities and 
coupling constants of strange particles directly 
from angular distributions is described. A rea- 
sonably detailed discussion is given of what can 
be gleaned by this method from the processes 
of associated production by photons on nucleons 
and by pions on protons, and from the absorp- 
tion and charge exchange scattering of charged 
K-mesons on nucleons. If the relative parity of 
K* to K° is odd then strange particle parities 
may be obtained from all these processes, while 
if this relative parity is even only photoproduc- 
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tion processes may be used to determine parities. 
Once these relative parities have been determined 
then coupling constants may be obtained from all 
the processes. 


HYPERFINE SPLITTING EFFECTS IN THE CAP- 
TURE OF POLARIZED »~ MESONS. H. therall, 
Carnegie Institute of Technology, Pittsburgh, 
Pennsylvania (Received January 16, 1959). 


We investigate the effect of the hyperfine split- 
ting of y-mesonic atom ground states on the neu- 
tron asymmetry from muon capture in hydrogen, 
in deuterium, and in complex nuclei with spin, 
which are treated by a one-particle model. It is 
shown that this can provide more information on 
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the capture interaction than the neutron asymmetry ayy) 


from spinless nuclei. Muon polarizations and 
gyromagnetic ratios in the hyperfine states are 
also discussed. 


FIXED-ANGLE DISPERSION RELATIONS FOR 
NUCLEON-NUCLEON SCATTERING. M. Cini, 
Istituto di Fisica e Scuola di Perfezionamento in 
Fisica Nucleare dell’ Universita, Roma, Italia, 
and Istituto Nazionale di Fisica Nucleare, Sezione 
di Roma, Roma, Italia, and S. Fubini and A. 
Stanghellini, CERN, Geneva, Switzerland (Re- 
ceived January 16, 1959). 


Starting from the same spectral representa- 
tion of the nucleon-nucleon scattering amplitude 
used by Chew in his determination of the pion- 
nucleon coupling constant, dispersion relations 
in which the scattering angle @ is kept constant 
are derived. An extrapolation procedure based 
on these dispersion relations is proposed with 
the purpose of obtaining an independent deter - 
mination of the pion-nucleon coupling constant 
and the possibility of discriminating between the 
different available sets of phase shifts. 

A first determination of the pion-nucleon coup- 
ling constant from these fixed-angle dispersion 
relations is attempted by means of a simplified 
extrapolation procedure, closely analogous to the 
effective-range approximation, and a result in 
good agreement with the currently accepted value 
is obtained. 
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